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1.  INTRODUCTION  AND  RESEARCH  OBJECTIVES  . 

...  /''•(  '  _ 

-  This  report  describes  work  carried  out  at  the  Gas  Turbine  Laboratory  at 

A 

MIT,  as  part  of  our  multi-investigator  effort  &n  high  performance  turbomachin- 
ery  fluid  dynamics.  Support  for  thif  program  is  provided  by  the  Air  Force 
Office  of  Scientific  Research  under  Contract  Number  F49620-85-C-0018,  Dr.  J.D. 
Wilson,  Program  Manager. 

The  present  report  gives  a  short  summary  of  the  work  for  the  period 
10/19/85  -  10/18/86.  For  further  details  and  background,  the  referenced  labora¬ 
tory  reports,  publications,  the  previous  final  report  (Reference  1)  covering 
the  period  10/1/81  -  9/30/84,  and  the  two  recent  progress  reports  on  the 
current  contract  should  be  consulted. 


Within  the  general  topic,  four  separate  tasks-  are  specified, 
in  brief: 


These  are. 


I.  Loss  mechanisms  and  loss  migration  in  transonic  compressors,  including 
development  of  advanced  instrumentation  for  measurements  of  wake  radial 
transport  and  analysis  of  unsteady  vortical  wake  structures.  £ 

II.  Experimental  and  theoretical  study  of  flows  in  casing  and  hub  treatment, 
including  mechanisms  for  stability  enhancement  in  compressors  and 
unsteady  fluid  dynamic  interactions  between  passage  and  groove  flows. 

III.  Computational  techniques  for  turbomachinery,  including  inverse  (design) 
calculation  procedures  for  transonic  turbomachine  blades  accounting  for 
viscid/ invisc id  interaction.  ' 

IV.  Theoretical  modelling  of  stability  and  unsteadiness  in  transonic  compres¬ 
sor  flow  fields,  including  analyses  of  unsteady  temperature  fluctuations 
due  to  vortex  shedding. 

The  work  carried  out  in  each  of  the  tasks  will  be  described  in  the  next  section. 


^  VW  »  *  .Ca  uV  —  V  aS.  a  a  a*»  a\  a  "a  — *a  a*a  CU  ’J  *  -  *  ."L-  l!  *  *a  .  *.  *.  I 
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TASK  I:  LOSS  MECHANISMS  AND  LOSS  MIGRATION  IN  TRANSONIC  COMPRESSORS 
Objectives 

The  primary  objective  of  this  effort  is  to  elucidate  and  quantify  the 
physical  mechanisms  which  produce  loss  (entropy)  in  a  high  speed  compressor, 
with  particular  attention  to  three-dimensional  and  unsteady  effects.  The 
motivations  for  this  work  are  twofold.  The  first  is  the  desire  to  dramati¬ 
cally  improve  the  design  point  performance  of  high  speed  compressors  to 
meet  the  ambitious  Air  Force  propulsion  goals  of  the  1990 's.  The  second  is 
the  realization  (spurred  by  AFAPL  data  showing  apparent  efficiencies  above 
100%)  that  there  are  important  physical  phenomena  influencing  the  behavior 
of  compressors  which  are  not  at  all  understood  or  accounted  for  in  present 
design  methodologies. 

Progress  During  This  Year 

Effort  during  the  past  year  has  been  concentrated  in  three  areas: 
analysis  of  the  effects  of  wake  structure  on  compressor  efficiency  and  the 
interpretation  of  performance  measurements;  an  ab  initio  calculation  of  the 
unsteady  flow  in  transonic  blade  rows  using  a  2-D,  time  accurate,  Reynolds 
averaged  Navier-Stokes  procedure;  and  continued  work  in  the  mapping  of  the 
radial  flow  in  a  transonic  compressor  rotor. 

Wake  Structure  Effects 

The  analysis  of  the  effects  of  wake  structure  on  apparent  compressor 
efficiency  is  a  follow-up  to  work  reported  earlier  which  identified  vortex 
street  structures  in  the  wakes  of  high  performance  transonic  compressors, 
through  both  experimental  and  analytical  investigations.  It  has  been  shown 
that  a  vortex  street  structure  in  the  compressor  wake  can  influence  the  com¬ 
pressor  efficiency,  noise,  and  structural  integrity.  The  analysis  also  shows 
that  the  wake  structure  can  influence  the  measurement  of  compressor  performance. 


introducing  artifacts  even  with  "perfect”  instrumentation.  For  the  stage 
studied  here,  this  measurement  error  is  of  the  order  of  1%  in  efficiency. 


Another  interesting  effect  of  the  vortex  street  is  to  redistribute  ther¬ 
mal  energy  in  the  compressor,  thus  introducing  total  temperature  gradients. 

The  total  temperature  gradients  combine  with  the  shock  structure  in  the 
compressor  to  yield  an  unusual  result,  adiabatic  efficiency  (the  common 
measure  of  compressor  performance)  and  entropy  change  are  no  longer  congruent 
as  they  are  in  flow  with  uniform  total  temperature.  Thus,  our  well  developed 
"feel"  for  these  flow  field  descriptors  is  erroneous.  Although  this  is  really 
a  bookkeeping  problem,  it  may  have  an  important  influence  in  how  we  design  and 
compare  turbomachines. 

Numerical  Flow  Field  Simulation 

To  further  investigate  compressor  wake  structure,  a  2-D,  time  accurate, 
viscous,  Reynolds  averaged  Navier-Stokes  calculation  of  a  rotor  passage  was 
performed  using  the  midspan  geometry  of  an  experimentally  investigated  compres 
sor.  This  ab  initio  calculation,  which  has  a  uniform  inflow  condition,  shows 
vortex  shedding  very  similar  to  that  observed  experimentally.  The  mean  fre¬ 
quencies  of  shedding  from  the  calculation  and  those  inferred  from  two  indepen¬ 
dent  measurement  techniques  are  ail  about  15  kHz,  3  times  the  rotor  passing 
frequency,  a  very  strong,  mutual  confirmation. 

A  surprising  result  of  the  calculation  was  a  low  frequency  (350  Hz) 
modulation  of  vortex  shedding,  apparently  caused  by  motion  of  the  separation 
point  along  the  suction  surface  of  the  airfoil.  This  oscillation  appears 
similar  in  nature  to  a  diffuser  instability  and  was  observed  (but  remained 
unexplained)  in  experimental  measurements.  The  importance  of  this  low  fre¬ 
quency  instability  is  not  fully  understood  but  since  a)  the  resultant  total 
pressure  fluctuations  are  on  the  order  of  5%  of  the  stage  pressure  rise,  and 


%  .  **. 


b)  the  instability  results  in  a  30%  fluctuation  in  blade  moment  at  a  frequency 


within  5%  of  the  first  blade  bending  frequency,  we  believe  that  further  inves¬ 
tigation  is  warranted. 

Overall,  from  this  work  we  have  made  the  following  observations. 

1.  The  wakes  of  high  Reynolds  number,  transonic  compressor  blades  can  consist 
of  shed  vortex  streets.  This  is  confirmed  by  measurement,  modelling,  and 
numerical  simulation. 

2.  The  vortex  shedding  scales  with  wake  thickness. 

3.  The  shedding  frequency  and  strength  are  sensitive  to  the  environment. 

4.  The  vortex  street  can  depress  the  average  wake  temperature. 

5.  The  effect  of  the  wake  structure  on  temperature  and  pressure  levels 
increases  with  the  square  of  the  freestream  Mach  number. 

6.  Wake  structure  can  influence  the  measurement  of  compressor  efficiency. 

7.  Changes  in  wake  structure  can  be  mistaken  for  changes  in  compressor 
performance . 

8.  The  wake  structure  can  drive  the  blade  shock  system,  inducing  loss. 

9.  The  wake  structure  can  be  influenced  by  other  fluid  dynamic  instabilities 
present  in  a  blade  passage  which  themselves  may  be  important  for  perform¬ 
ance  and  structural  dynamic  reasons. 

10.  Caution  is  required  when  interpreting  steady  state  solutions  of  inherently 
unsteady  flowfields. 

More  detail  is  included  in  the  paper  preprint  appended. 

Radial  Flow  Migration  in  a  Transonic  Rotor 

The  structure  of  the  radial  flows  in  a  turbomachine  can  have  an  important 
influence  on  the  behavior  of  a  given  stage  as  well  as  on  the  performance  of 
subsequent  blade  rows  downstream.  These  radial  flows  give  rise  to  spanwise 


mixing  which  is  a  key  determinant  of  multistage  compressor  performance.  The 
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nature  of  magnitude  of  this  process  has  recently  become  a  topic  of  considerable 
debate  (Atkins  and  Smith  versus  Gallimore  and  Cumpsty)  fueled  by  its  very  real 
practical  importance.  In  all  cases  to  date,  these  radial  mixing  processes 
have  been  treated  only  in  the  time  average.  In  this  work,  we  are  taking  time 
resolved  measurements  of  the  flow  field  using  a  tracer  gas  to  elucidate  the 
radial  flow  transport  in  the  rotor.  The  motivation  is  both  to  unsort  the 
radial  flow  structure  (spanwise  vortex  coherence,  flows  in  boundary  layer 
seapration  regions,  etc.),  and  to  quantify  the  radial  mixing  mechanisms. 

The  mapping  of  the  radial  flow  in  the  AFAPL  high  thru-flow  rotor  is  being 
carried  out  using  gas  injection  in  the  MIT  Blowdown  Compressor  Facility. 
Basically,  a  circumferential  sheet  of  gas  (helium)  is  injected  upstream  of  the 
rotor  and  the  time  resolved  concentration  distribution  of  the  injectant  is 
measured  at  the  rotor  outlet.  The  experimental  technique  has  been  developed 
during  previous  years  of  this  effort.  The  layout  of  the  apparatus  is  shown  in 
Figs.  1  and  2.  The  injector  was  first  run  in  the  tunnel  without  the  compressor 
stage  to  establish  a  baseline  mixing  rate  for  the  injectant  sheet.  The  stage 
was  then  run  with  the  injector  present,  with  and  without  the  injectant  turned 
on,  to  indicate  the  effect  of  the  injector  hardware  (blockage  and  viscous 
wake)  on  the  stage  performance.  The  experiment  with  injectant  was  then  run, 
sample  data  from  which  is  shown  in  Fig.  3. 

The  experiment  has  proven  considerably  more  difficult  than  first  expected, 
with  detection  sensitivity  the  principal  problem.  For  this  reason,  the  t -acer 
gas  has  been  changed  from  CO2  (the  best  match  to  the  main  flow)  to  helium  (ten 
times  the  sensitivity).  To  date,  the  flow  has  been  mapped  with  injection  at 
the  midspan  and  hub  positions. 

It  is  planned  to  finish  the  experimental  mapping  over  the  course  of  the 
next  year  and  to  reduce  and  model  the  data.  Preliminary  analysis  shows  the 
time  averaged  results  match  those  of  Gallimore  and  Cumpsty. 


Figure  1:  Scale  drawing  of  the  test  section 


figure  .! :  Top  view  of  the  experimental  arrangement 
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VORTEX  SHEDDING  IN  HIGH  SPEED  COMPRESSOR  BLADE  HAKES 

A.H.  Epstein,  J.B.  Gertz,  P.R.  Owen  and  M.B.  Giles 
Gas  Turbine  Laboratory 
Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139  USA 

ABSTRACT 

The  wakes  of  highly  loaded  compressor  blades  are  generally  considered  to 
be  turbulent  flows.  Recent  work  has  suggested  that  the  blade  wakes  are  domi¬ 
nated  by  a  vortex  street-like  structure.  The  experimental  evidence  supporting 
the  wake  vortex  structure  is  reviewed.  This  structure  is  shown  to  redistribute 
thermal  energy  within  the  flowfield.  The  effect  of  wake  structure  on  conven¬ 
tional  aerodynamic  measurements  of  compressor  performance  is  noted.  A  two- 
dimensional,  time  accurate,  viscous  numerical  simulation  of  the  flow  exhibits 
both  vortex  shedding  in  the  wake  and  a  lower  frequency  flow  instability  which 
modulates  the  shedding.  The  numerical  results  are  shown  to  agree  quite  well 
with  the  measurements  from  transonic  compressor  rotors. 

1 .  INTRODUCTION 

Although  bluff  bodies  and  low  speed  airfoils  have  long  been  known  to  shed 
vortices  into  their  wakes,  the  wakes  of  high  speed,  high  Reynolds  number 
turbomachinery  blading  have  generally  been  considered  to  be  turbulent  and 
unstructured.  Usually,  the  wake  is  described  only  in  terms  of  a  time  averaged 
velocity  profile.  This  is  consistent  with  the  view  that,  assuming  uniform 
inflow  and  excepting  turbulence,  the  flow  in  the  frame  moving  with  the 
compressor  rotor  is  uniform  and  that  variations  observed  in  the  stationary 
frame  are  primarily  due  to  blade  to  blade  geometric  differences.  This  view, 
however,  is  erroneous  in  detail  with  practical  import  for  compressor  design. 

There  is  a  considerable  body  of  experimental  data  which  shows  that  blunt 
trailing  edge  airfoils  typical  to  turbines  shed  vortex  streets  [1],  [2],  [3]. 


Recent  work  has  shown  that  a  similar  phenomenon  may  occur  in  sharp  trailing 
edge  transonic  compressors  as  well  [4]. 

If  we  assume  for  the  moment  that  high  speed  compressor  blade  wakes  can 
consist  of  a  vortex  street-like  structure,  then  a  number  of  questions 
immediately  become  relevant.  These  include:  do  all  blade  wakes  contain 
vortices;  why  are  the  vortices  rarely  observed  experimentally;  how  are  the 
vortex  streets  formed  and  what  is  their  structure;  and  what  is  the  practical 
importance  of  the  wake  structure  to  the  compressor  designer?  We  will  address 
the  last  question  first  so  as  to  provide  a  groundwork  for  the  discussions 
which  follow. 

In  the  context  of  high  efficiency,  high  performance  transonic  compressors 
and  fans,  the  influence  of  the  wake  temporal  and  spatial  structure  for  a  given 
time  averaged  velocity  defect  can  be  surprisingly  large.  The  wake  structure 
can  influence  the  compressor  aerodynamics,  noise,  and  structural  integrity. 
Vortex  shedding  decreases  the  blade  row  efficiency  (higher  base  drag  anc  shock 
loss,  increased  unsteadiness)  and  induces  artifacts  m  standard  measurement 
techniques.  An  increased  level  of  unsteadiness  into  following  Diade  rows  can 
increase  the  loss  there  and  alter  the  mass  flow  characteristics,  independently 
of  whether  that  blade  row  is  itself  shedding  vortices.  The  interaction  of  the 
shed  vortex  street  with  subsequent  blade  rows  is  a*sc  a  source  of  noise. 
Obviously,  the  unsteady  loading  due  to  tne  vortex  streets  and  related  phenomena 
can  have  deleterious  effects  on  the  structura*  integrity  of  a  compressor  when 
the  aerodynamic  excitation  frequencies  coincide  witr.  tnose  which  are  structur¬ 
ally  important.  The  blade  waxe  structure  a~so  influences  compressor  design 
and  analysis  in  that  it  represents  a  physica*  pnenomenon  not  usually  modelled 
(either  because  the  analysis  is  inviscid  or  steady  state).  Thus,  the  design 
intent  may  not  represent  a  physically  realizable  system. 


This  paper  is  an  extension  of  previous  work  in  which  the  presence  of 
vortex  streets  was  xnferred  from  experimental  measurements  and  their  effect  on 
the  blade  shock  system  discussed  [4],  [5].  Here,  we  review  that  work,  discuss 
the  effects  of  vortex  streets  on  the  temperature  distribution  within  the  flow, 
note  the  influence  of  the  wake  structure  on  conventional  aerodynamic  perfor¬ 
mance  measurements,  present  the  results  of  an  ab  initio  numerical  simulation 
showing  similar  behavior,  and  suggest  areas  for  further  investigation. 

2.  REVIEW  OF  EXPERIMENTAL  OBSERVATIONS 

Flow  visualization  studies  have  demonstrated  the  existence  of  vortex 
streets  in  the  wakes  of  transonic  flat  plate  and  turbine  airfoil  cascades  [1], 
[2],  [3]  but,  excepting  the  special  case  of  acoustic  resonance  enhanced 
shedding  [6],  [7],  shedding  in  high  speed  compressor  airfoils  had  not  been 
conclusively  reported.  Hot  wire  measurements  in  the  wakes  of  controlled 
diffusion  airfoil  cascades  yielded  ambiguous  results  [8]  with  one  configuration 
showing  evidence  of  periodicity  in  the  wake  while  another  did  not.  In  the 
context  of  hydrofoils,  Blake  [9]  suggested  that  all  bodies  shed  but  may  do  so 
in  a  discontinuous  fashion,  i.e.,  shed  in  bursts.  Lack  of  phase  coherence 
between  the  bursts  would  then  tend  to  obscure  the  shedding  when  examined  with 
spectral  analysis  techniques.  This  is  consistent  with  the  observation  tMt 
the  least  ambiguous  information  has  been  instantaneous  flow  visualizations. 

Time  resolved  measurements  in  the  outflow  of  a  transonic  compressor 
rotor  operating  near  the  compressor's  peak  efficiency  points  revealed  high 
frequency  ( 3  to  4  times  blade  passing)  total  pressure  and  temperature 
fluctuations  of  substantial  amplitude  in  the  core  flow  between  the  wakes,  as 
well  as  large  fluctuations  in  the  wake  strength  of  any  particular  blade  [5], 
Fig.  1.  Similar  observations  were  made  in  three  different  single  stage 
transonic  compressors  in  three  different  test  facilities  using  a  variety  of 


high  frequency  response  probe  types  [11].  The  machines  differed  in  design 
intent  (both  commercial  fan  and  high  speed  military  designs  were  included), 
but  all  operated  at  high  rotor  efficiency  (above  90%).  The  outflow  of  all 
three  rotors  demonstrated  similar  types  and  degrees  of  rotor  relative 
unsteadiness,  the  unsteadiness  being  maximum  near  the  design  or  maximum 
efficiency  operating  point.  This  indicates  that  rotor  relative  unsteadiness 
is  a  common  phenomenon  in  transonic  compressors. 

A  laser  anemometer  physically  measures  one  or  more  components  of  the 
instantaneous  velocity  at  a  point  in  space.  Because  of  the  relatively  high 
level  of  turbulence  common  to  turbomachinery  flowfields,  the  instantaneous 
velocity  measurements  are  usually  averaged  and  the  result  presented  as  the 
average  velocity  at  that  point  (Fig.  2a).  The  velocity  statistics  at  the 
point  can  also  be  examined  however  (although  this  is  not  commonly  done),  and 
can  be  presented  in  terms  of  a  probability  density  distribution  (PDD),  a 
histogram  showing  the  number  of  observations  made  at  each  value  of  velocity. 
For  a  turbulent  flow,  the  shape  of  this  histogram  should  be  Gaussian.  Laser 
anemometer  measurements  in  the  vicinity  of  the  passage  shock  have  shown  a 
bimodal  velocity  distribution  indicative  of  a  shock  moving  to  either  side  of 
the  measurement  point  [10].  These  observations  were  explained  in  terms  of 
small  (0.5%  of  axial  chord)  motion  of  the  passage  shock  about  its  mean 
position.  Vortex  shedding  in  the  blade  wake  was  advanced  as  a  possible  driver 
of  this  motion  [5], 

Subsequent  laser  anemometer  measurements  in  the  rotor  wakes  revealed  that 
the  velocity  distributions  were  bimodal  there  as  well;  i.e.,  two  velocities 
were  equally  likely  with  few  measurements  observed  at  the  "average"  velocity 
as  it  would  be  normally  derived  from  such  data  [4],  Fig.  2.  Simultaneous  time 


resolved  temperature  and  pressure  probe  measurements  in  the  rotor  outflow 
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snowed  unsteadiness  relative  to  the  rotor  at  two  time  scales.  One  was  the 
high  frequency  disturbances  mentioned  above.  The  second  appeared  primarily  as 
a  modulation  of  the  wake  flow  with  frequency  components  on  the  order  of  1/2  to 
3  times  shaft  rotational  speed.  This  modulation  could  not  be  explained  in 
terms  of  blade  to  blade  geometric  differences  since  the  fluctuations  were 
primarily  aperiodic  with  rotor  rotation  (and  thus  not  locked  to  the  geometry) 
and  were  several  times  the  magnitude  of  the  periodic  disturbances.  A  striking 
feature  of  the  disturbances  themselves  was  the  large  total  temperature  fluctu¬ 
ations  in  the  wakes — ±3-5%  of  the  mean  total  temperature  (i.e.,  10-20°C)-- 
indicating  local  regions  of  intense  cooling  and  heating.  The  instantaneous 
adiabatic  efficiency  calculated  from  these  measurements  showed  a  concomitant 
variation  with  some  local  regions  appearing  as  over  100%.  These  temporal 
fluctuations  and  laser  anemometer  measurements  were  explained  as  resulting 
from  a  vortex  street  structure  in  the  blade  wakes.  The  over  100%  efficiency 
observations  have  yet  to  be  explained. 

3.  DESCRIPTIONS  OF  THE  BLADE  WAKE  AND  VORTEX  STREET 

In  the  previous  work,  the  wake  was  modelled  as  two  staggered  rectilinear 
rows  of  Rankine  vortices  of  opposite  sign  in  a  uniform  freestream  [11).  The 
vortices  consist  of  an  inner  region  with  a  forced-vortex  core  and  an  outer 
region  following  the  irrotational  flowfield  of  a  classic  von  Karman  vortex 
•street.  This  model  was  fit  to  the  laser  anemometer  data  so  that  the  time 
averaged  velocity  profile  and  the  statistical  distribution  of  velocities  (the 
probability  density  distribution)  in  the  wake  would  match.  The  vortex  size 
and  strength  were  adjusted  to  fit  the  velocity  profile  while  the  ratio  of  tne 
streamwise  vortex  spacing  (A)  to  the  distance  between  vortex  rows  (H)  deter¬ 
mined  the  velocity  statistics  (see  Fig.  5). 

As  can  be  seen  in  Fig.  3,  the  model  fits  the  data  relatively  well,  which 
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is  a  consistency  check.  The  model  prediction  of  a  shedding  frequency  of  16+2 
kHz  is  quite  close  to  the  14-15  kHz  inferred  from  the  core  flow  pressure 
fluctuations.  Furthermore,  the  model  readily  explains  the  high  level  of 
fluctuations  observed  in  the  wake  with  the  high  frequency  response  probes. 

This  is  an  artifact  of  sampling  caused  by  the  random  position  of  the  vortices 
in  the  wakes  as  the  compressor  revolves  past  the  probe  location.  Figure  4 
compares  the  absolute  frame  rotor  outflow  total  pressure  fluctuations  predicted 
by  the  model  with  measurements. 

While  the  vortex  street  model  in  [11]  does  a  good  job  of  explaining  many 
of  the  experimental  observations,  it  is  a  static  model,  describing  the  wake 
state  only  at  a  particular  axial  station.  Thus,  it  contains  no  information  on 
the  vortex  formation  process,  its  evolution  or  decay.  At  this  time  we  know  of 
no  published  analytical  model  describing  this  process  in  detail. 

4.  DYNAMIC  ENERGY  REDISTRIBUTION 

One  of  the  more  interesting  implications  of  identifying  the  structure  of 
the  wake  of  a  transonic  compressor  blade  as  containing  a  vortex  street  is  that 
the  unsteady  pressure  field  in  the  blade  relative  frame  can  redistribute 
thermal  energy.  The  vortex  street  propagates  at  a  velocity  different  from 
that  of  the  mean  flow  and  entrains  fluid  into  and  through  the  wake.  As  the 
entrained  fluid  passes  through  the  pressure  field  of  the  vortex,  its  total 
enthalpy  can  oe  changed.  In  essence,  the  vortex  street  can  be  thought  of  as 
an  array  of  tiny  turbomacr.mes : 

Tne  dynamics  of  tr.is  proo^err  nave  oeen  addressed  ov  McCune  [11  ]  witr.  a 
straightf orwarc  unsteady  energy  equation  analysis.  He  snows  tnat  tne  pressure 
f -eld  fluctuations  from  tne  vortices  cnanae  tne  tota-  temperature  distribution 
in  tne  fluid  anc  tnat  tne  amplitude  of  this  effect  scai.es  with  tne  square  of 
tne  freestream  Macr.  number.  At  Macr.  1,  tne  vortices  are  predicted  to  produce 
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a  temperature  fluctuation  of  3-5%  of  the  freestream  total  temperature,  quite 
close  to  the  observed  value.  This  scaling  with  Mach  number  explains  why  the 
phenomenon  was  first  noticed  in  a  high  speed  machine  (tip  relative  Mach  number 
of  2.2). 

Kurosaka  et  al.  treat  the  problem  in  terms  of  the  general  Eckert-Weise 
effect  [13]  in  which  the  recovery  factor  at  the  rear  of  a  shedding  right 
circular  cylinder  can  be  negative  and  the  wake  centerline  cooled  to  below  the 
freestream  inflow  value,  showing  that  this  is  the  result  of  vortex  street 
behavior  [14].  A  key  point  is  that  this  cooling  by  the  vortices  is  observed 
with  steady  state  instrumentation,  i.e.,  the  average  centerline  temperature  is 
depressed.  There  must,  of  course,  be  a  concomitant  increase  in  total  tempera¬ 
ture  in  the  surrounding  flow  outside  the  wake.  Kurosaka  also  demonstrated 
that  the  vortex  shedding  and  temperature  separation  can  be  significantly 
enhanced  with  acoustic  feedback,  essentially  extending  acoustic  resonance  work 
to  include  energy  separation  [6],  [7], 

The  intense,  localized  hot  and  cold  spots  generated  by  the  vortex  street 
explain  the  large,  abrupt  temperature  fluctuations  observed  in  the  time 
resolved  total  temperature  measurements  of  the  compressor  wakes. 

5.  WAKE  STRUCTURE  EFFECTS  ON  APPARENT  COMPRESSOR  EFFICIENCY 

An  earlier  work  on  high  frequency  passage  shock  motion  driven  by  vortex 
shedding  identified  three  mechanisms  by  which  the  loss  in  the  compressor  could 
increase  [5].  The  first  is  the  increase  in  entropy  rise  across  a  shock  wave 
undergoing  small  periodic  axial  motion  compared  to  that  for  a  steady  shock  at 
the  same  average  approach  Mach  number.  This  change  was  quite  small,  amounting 
to  only  a  0.1  to  0.2  percent  decrease  of  adiabatic  efficiency  in  the  rotor 
studied.  The  second  loss  source  was  the  small  scale  nonuniformities  in  total 
pressure  generated  by  the  oscillating  shock  wave.  In  this  case,  it  was 
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assumed  that  due  to  the  small  spatial  extent  of  the  perturoation  (l/8^1/4 
chord),  the  enthalDy  would  not  be  recovered  as  pressure  rise  in  the  diffusion 
process  in  the  stator  but  rather  appear  as  a  mixing  loss  downstream.  This 
loss  was  calculated  to  be  of  the  same  order  as  the  wake  mixing  loss,  about  1% 
(10%  of  the  total  measured  stage  loss).  The  third  loss  mechanism  is  the 
amplification  of  the  spatial  nonuniformities  in  the  stage  outflow,  discussed 
above,  by  the  shock  system  in  a  following  transonic  stage.  This  loss  was 
estimated  to  be  2  to  3  times  as  great  as  that  in  the  first  stage,  i.e.,  2  to  3 
percent  of  stage  adiabatic  efficiency. 

In  this  section  we  wish  to  discuss  not  additional  thermodynamic  loss 
mechanisms  (i.e.,  production  of  entropy),  but  rather  apparent  losses,  artifacts 
induced  in  the  measurement  processes  and  their  interpretation  by  the  periodic 
nature  of  the  wake  structure.  If  not  properly  accounted  for,  these  artifacts 
can  result  in  an  erroneous  estimate  of  compressor  performance. 

Before  doing  so,  however,  we  wish  to  point  out  a  simple  curiosity.  Adia¬ 
batic  efficiency  is  often  used  to  express  turbomachine  loss  (entropy  produc¬ 
tion).  In  an  unsteady,  transonic  flow,  however,  the  change  in  adiabatic 
efficiency  need  not  be  congruent  with  the  entropy  production,  as  it  must  in  a 
steady  flow. 

To  illustrate  this,  we  will  define  the  local  adiabatic  efficency,  n ,  in 
the  usual  manner  relating  the  total  pressure,  Pt,  and  total  temperature,  Tt, 
ratios  in  the  absolute  ( laboratory )  frame  as  follows. 
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wnere  the  subscripts  I  and  I  denote  stations  upstream  and  downstream  of  the 
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rotor  olaae.  Station  1  is  assumed  to  have  uniform  conditions  of  total  tempe¬ 
rature  and  pressure  as  would  be  the  case  with  the  first  stage  of  a  compressor 
Station  2  consists  of  a  freestream  region  and  a  wake  region,  with  the  blade 
wakes  represented  by  vortex  streets  as  modelled  in  [11].  The  modelled  wake 
total  pressure  and  temperature  contours  are  shown  in  Figs.  5  and  6.  The 
change  in  entropy,  S,  from  station  1  to  station  2  may  be  expressed  as, 


Y-l 


VV/CP 

fV 

K 

Y 

abs 

'rn  ' 

t2 

T 

l  tlj 

abs 

This  can  be  seen  in  Fig.  7  for  a  flow  with  the  compressor  blade  wake  vortex 
model.  The  freestream  stagnation  temperature  ratio  is  1.175.  As  expected, 
the  entropy  variation  is  zero  in  the  regions  outside  the  vortex  cores  (since 
the  model  assumes  that  all  the  entropy  change  in  the  wakes  is  in  the  vortex 
cores  themselves). 

Equation  (1)  for  the  efficiency  may  now  be  re-written  by  eliminating  the 
stagnation  pressure  ratio  using  Eq.  (2). 
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Equation  (3)  shows  that  in  a  flow  of  uniform  total  temperature,  the  adiabatic 
efficiency  is  a  unique  function  of  the  entropy  change.  Conversely,  if  there 
is  no  change  in  entropy  from  station  1  to  station  2  then  the  efficiency  is 
exactly  1.0  regardless  of  the  variations  in  total  temperature.  However,  if 


between  stations  1  and  2  there  is  some  loss  mechanism  (such  as  a  normal  shock 
wave)  so  that  S2  >  in,  the  freestream  region,  then  the  efficiency  will  be 
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less  than  1.0  and,  in  fact,  will  vary  through  the  flow  as  the  total  temperature 
varies.  In  the  case  of  a  transonic  compressor  with  a  vortex  street  wake,  the 
efficiency  will  vary  even  outside  the  vortex  cores  in  the  wake  (which  contain 
all  the  entropy  change)  where  the  flow  is  otherwise  modelled  as  irrotational 
and  inviscid,  i.e.,  the  value  of  (S2“S})  is  constant  (as  shown  in  Fig.  7). 

Since  the  flow  in  the  absolute  frame  is  unsteady  (with  or  without  vortex  shed¬ 
ding)  and  the  gradient  of  entropy  is  zero  outside  the  cores,  the  total  tempera¬ 
ture  must  vary  according  to  Crocco's  theorem  which  may  be  written  as  follows, 

— ► 

TVS  =  V(c  Tt)  -  q(Vxq)  +  |3  (4) 

In  the  regions  outside  the  vortex  cores  this  reduces  to. 
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Equation  (5)  implies  that  the  total  pressure  must  vary  as  the  total  temperature 
changes  so  as  to  keep  the  entropy  (S2~S]_)  constant.  Iq  any  case,  since  (S2-S}) 
is  constant  and  T^/Tti  varies  throughout  the  vortex  street,  the  efficiency 
will  vary  as  well,  even  though  the  entropy  is  constant!  This  can  be  seen  by 
comparing  Figs.  7  and  8.  The  efficiency  is  not  congruent  with  the  entropy 
change.  Because  of  the  energy  separation  mechanisms  discussed  in  Section  4, 
the  total  temperature  in  local  regions  can  be  greater  or  less  than  that  of  the 
freestream.  Thus,  there  are  regions  in  which  the  efficiency  is  greater  than 
the  freestream  value.  It  should  be  pointed  out  with  reference  to  Eq.  (3) 
that,  so  long  as  the  entropy  increases,  a  variation  in  total  temperature  cannot 
result  in  an  adiabatic  efficiency  above  one. 

It  is  important  to  note  that  this  difference  between  entropy  rise  ar.c 
efficiency  requires  only  that  the  flow  have  a  change  in  entropy  Between  t.ne 
measurement  stations  (e.g.,  a  shock  wave)  and  a  spatial  variation  ir.  total 
temperature.  Unsteadiness  plays  a  part  only  in  that,  in  the  absence  of  neat 


conduction  anc  viscous  stresses  ,  1j  ]  ,  a  variation  ir.  total  temperature  can 
only  be  produced  by  unsteady  pressure  forces.  The  unsteadiness  is  required  in 
the  absolute  frame,  not  the  blade  relative  frame,  and  thus  a  blade  relative 
unsteady  wake  is  not  a  prerequisite.  Rather,  this  effect  is  common  to  all 
transonic  and  supersonic  turbomachines .  Other  than  suggesting  that  adiabatic 
efficiency  is  a  misleading  measure  of  turbomachine  performance,  the  practical 
significance  of  this  observation  has  yet  to  be  enunciated,  but  the  authors 
find  it  interesting. 

The  structure  of  the  compressor  blade  wakes  can  make  a  difference  in  the 
aerodynamic  efficiency  as  commonly  measured  with  aerodynamic  probes  in  the 
rotor  outflow.  To  illustrate  this,  we  will  consider  two  transonic  compressor 
rotors  with  the  same  geometry,  true  mass  averaged  total  pressure  and  tempera¬ 
ture  rise,  and,  therefore,  the  same  mixed  out  adiabatic  efficiency.  One  blade 
will  be  assumed  to  have  a  classic  turbulent  wake  while  the  other's  wake  will 
consist  of  a  vortex  street  with  the  same  average  velocity  defect.  (Here  we 
use  the  term  turbulent  rather  loosely,  meaning  random,  thus  ignoring  any  turbu¬ 
lent  structure  which  may  exist.)  In  other  words,  the  wakes  as  measured  by 
conventional  laser  anemometry  techniques  would  appear  identical.  Conventional, 
low  frequency  response  prooes  are  then  used  to  accurately  measure  the  true 
time  average  of  the  total  temperature  and  pressure  at  the  rotor  outflow.  (Since 
this  is  a  two-dimensional  analysis,  the  area  averaged  quantities  measured  by 
prooes  are  actually  line  averages.)  The  probes  are  assumed  to  De  free  from 
dynamic  effects.  We  will  taxe  the  average  flow  conditions  at  the  rotor  blade 
row  exit  plane  to  oe  tnose  measured  ir.  a  transonic  rotor  [4],  Table  1.  Note 
that  the  total  temperature  _s  r.ig.ner  ir.  the  flow  outside  the  wakes  ir.  the 
vortex  street  case  compared  to  tnat  ir.  tne  turbulent  wane  case.  This  is 
required  since  the  vortex  street  cepresses  the  wane  centerline  temperature  but 
the  total  energy  flux  out  of  tne  D.aae  row  must  oe  tne  same  ir.  Doth  cases. 


TABLE  1 

ROTOR  EXIT  CONDITIONS  ASSUMED  FOR  MIXING  CALCULATION 


Freestream 


Vortex  Wake 


Turbulent  Wake 


Total  pressure  ratio 
Total  temperature  ratio 
Absolute  exit  flow  angle 
Wake  width/Passage  width 


1.646 

1.175 

45° 

0.260 


1.646 

1.170 

45° 

0.260 


Table  2  shows  the  resultant  average  flow  condition  calculated  with  the 
two  wake  states.  Results  are  given  both  close  to  the  rotor  (the  measurements 
were  made  at  130%  of  the  axial  chord)  where  the  wake  structures  are  distinct, 
and  far  downstream  where  the  wakes  have  totally  mixed  out.  The  mass  averaged 
quantities  agree  at  both  stations  of  course.  The  time  averages  close  to  the 
rotor  do  not,  however.  The  time  average  of  the  vortex  street  wake  case  is 
quite  close  to  the  mass  average  (probably  by  coincidence)  but  the  turbulent 
wake  case  efficiency  is  measured  at  over  1%  low.  (The  difference  between  time 
or  area  and  mass  averages  is  well  known.  The  point  is  that  the  low  frequency 
response  probes  can  only  measure  the  time  average.)  Far  downstream,  the  mass 
averaged  efficiency  has  decreased  by  1/2%  (which  is  why  probe  stations  are 
normally  placed  as  far  downstream  as  possible). 


TABLE  2 

ROTOR  OUTFLOW  AVERAGE  LABORATORY  FRAME  CONDITIONS 


Near 

Rotor  ( 130%  Axial  Chord) 

Vortex  Wake 

Turoulent  Wake 

Mass 

averaged  total  pressure  ratio 

1.644 

1.642 

Mass 

averaged  total  temperature  ratio 

1.175 

1 . 175 

Mass 

averaged  adiabatic  efficiency 

0.873 

0.871 

Prooe 

indicated  total  pressure  ratio 

1 . 644 

1.641 

ProDe 

indicated  total  temperature  ratio 

1  . 175 

Probe 

indicated  adiabatic  efficiency 

0.872 

0.859 

Far  Downstream  (Mixed  Out) 


Total  pressure  ratio 

1.648 

1.645 

Total  temperature  ratio 

1  .  177 

1  .  176 

Adiabatic  efficiency 

0.866 

0.867 

If  the  wake  structure  were  always  the  same,  the  difference  m  measured 
efficiencies  due  to  the  structure  could  be  mixed  in  with  empirical  prone 
calibration  factors  and  thus  be  calibrated  out  (in  theory).  If,  nowever,  a 
small  change  in  blade  design  or  turbomachine  operating  point  were  to  alter  th< 
wake  structure  alone,  there  wou-d  be  an  erroneous  change  in  measured  aerodyna¬ 
mic  efficiency,  an  error  of  over  1%  for  the  stage  studied  here.  In  other 
words,  a  change  in  the  waxe  structure  can  appear  as  a  change  in  stage  effi¬ 
ciency,  even  though  tne  momentum  and  thermal  energy  flux  through  tne  machine 
(the  pressure  and  temperature  rise)  have  not  changed.  This  error  would 
disappear  if  the  measurements  were  made  with  sufficient  time  resolution  to 
resolve  the  wake  structure. 

Thus,  wake  structure  can  introduce  artifacts  in  the  measurement  of 
compressor  performance,  introducing  apparent  changes  not  representative  of  th< 
state  of  the  fluid  exiting  the  machine.  Because  the  magnitude  of  the 
temperature  separation  in  the  vortex  street  scales  with  the  square  of  the  Mac: 
number,  this  effect  should  be  most  important  for  transonic  and  supersonic 
turbomachines.  The  argument  should  hold  true  for  turbine  flows  as  well,  out 
this  has  not  been  verified. 

6.  NUMERICAL  SIMULATION  OF  COMPRESSOR  BLADE  WAKE  SHEDDING 

To  further  investigate  the  compressor  waxe  structure,  tne  flow  tr.rough  a 
single  blade  passage  was  simulated  with  a  computational  fluid  dynamic  'CFD) 
technique.  A  two-dimensional,  time  accurate,  Reynolds  averaged,  explicit 
Navier- Stokes  calculation  was  done  with  a  relatively  fine  grid  and  very  small 
time  step  to  insure  good  spat-a*  and  temporal  resolution.  Tne  midspar  airfci) 
geometry  cf  tne  transonic  rotor  measured  in  [4]  was  -sec.  Tc  acccnocate 
supersonic  inflow  limitations  m  tne  coae ,  tne  inlet  relative  Macr.  number  was 
reduced  from  tne  1.1“  of  tne  measurement  condition  tc  1.92  for  tne  calcuj.atioi 


Inlet  total  conditions  were  adjusted  to  keep  the  Reynolds  number  based  on 
axial  chord  at  1.2  *  10^.  The  inlet  flow  at  the  upstream  boundary  is  specified 
as  uniform.  Thus,  there  is  no  external  excitation  or  periodicity,  only  that 
generated  by  the  flow  itself.  More  details  on  the  calculation  can  be  found 
in  [15]. 

The  primary  ooservation  to  be  made  from  this  ab  initio  calculation  is 
that  vortices  are  shed  into  the  blade  wake.  Figure  9  is  a  plot  of 
instantaneous  velocity  vectors  in  the  blade  relative  frame  near  the  trailing 
edge.  It  shows  the  vortex  structure  in  the  wake  as  well  as  a  separated  region 
on  the  suction  surface.  (The  vector  length  denotes  the  magnitude  of  the  local 
flow  velocity,  the  orientation  the  flow  direction.)  Only  two  vortices  are 
readily  seen  here  due  to  the  change  in  vortex  core  translational  velocity  as 
the  vortices  propagate  downstream.  This  has  the  effect  of  blurring  the  flow 
structure  in  any  frame  other  than  that  moving  with  the  cores.  A  similar 
problem  in  visualizing  the  vortices  was  encountered  in  the  modelling  effort 
[11].  In  both  cases,  the  ready  identification  of  the  vortex  structure  is 
strongly  dependent  on  the  frame  of  reference  chosen,  making  the  vortex  street 
fairly  elusive  and  hard  to  see. 

The  shedding  periodicity  shows  up  quite  clearly  in  the  calculated  trailing 
edge  static  pressure.  Fig.  10.  Also  of  interest  is  the  lower  frequency 
present  which  modulates  both  the  amplitude  and  the  frequency  of  the  vortex 
shedding.  This  is  quite  similar  to  low  frequency  modulation  observed  in  the 
experimental  measurements  [11].  In  the  simulation,  the  lower  frequency 
correlates  with  the  motion  of  the  separation  point  along  the  suction  surface 
and  with  axial  motion  of  the  passage  shock,  and  thus  may  be  a  shock  boundary 
layer  interaction.  Fig.  11.  There  is  also  a  30%  fluctuation  in  blade  moment 


at  this  frequency,  which  is  low  enough  (^350  Hz)  to  be  of  concern  to  the 


structural  designer.  The  exact  cause  of  this  low  frequency  movement  of  the 
separation  point  is  not  yet  clear  but  appears  quite  similar  to  instabilities 
observed  in  high  speed  diffusers. 

The  modulation  of  the  vortex  shedding  frequency  by  the  low  frequency 
oscillation  was  considerable,  a  factor  of  two.  The  strength  of  the  vortices 
varied  inversely  with  the  shedding  frequency.  The  shedding  frequency  range 
from  the  CFD  calculation  is  compared  with  those  inferred  from  the  laser 
anemometer  probability  density  distributions  and  from  the  core  flow-shock 
motion  fluctuations  in  Fig.  12.  These  three  estimates  are  completely 
independent  and  show  quite  good  agreement.  The  large  fluctuations  in  fre¬ 
quency  shown  in  the  numerical  simulation  tend  to  explain  the  difficulty 
encountered  in  extracting  a  single,  unambiguous  frequency  estimate  from  the 
experimental  measurements.  Fluctuations  of  this  magnitude  may  also  blur  the 
bimodal  anemometer  histograms. 

The  blade  relative  total  pressure  in  the  numerical  simulation  as  would  be 
measured  with  a  fixed  laboratory  frame  probe  as  the  rotor  passes  is  compared 
with  measurements  [11]  in  Fig.  13.  The  qualitative  agreement  is  excellent. 

The  calculation  clearly  captures  the  wake  modulation  evident  in  the  measurements 

The  high  frequency  jitter  of  the  passage  shock  at  the  shedding  frequency 
inferred  in  [5]  is  not  observed  in  the  calculation.  However,  since  the 
predicted  shock  motion  is  no  more  than  one  grid  cell  size  and  the  numerics 
spread  the  shock  over  five  grid  points,  this  is  not  surprising.  A  calculation 
with  a  much  finer  grid  size  would  be  required  in  order  tc  address  tr.is  proc.err 
properly . 

Overall,  the  numerical,  CFD  simulation  agrees  extremely  well  with  the 
experimental  observations  and  the  analytical  model.  The  one  area  ir.  which  the 
numerical  simulation  does  not  add  information  is  the  decay  of  the  waxe 


structure  as  it  is  convected  downstream  since  the  numerical  damping  overwhelms 
most  physical  dissipation  mechanisms. 

7.  DISCUSSION  AND  CONCLUSIONS 

Considerable  effort  has  traditionally  been  spent  on  establishing  the 
proper  parameter  with  which  to  correlate  vortex  shedding.  A  Strouhal  number 
based  on  trailing  edge  thickness  is  commonly  used,  especially  for  blunt 
trailing  edge  oodies  [2],  [3],  [9].  For  compressor  blades,  the  wake 
displacement  thickness  has  been  suggested  as  a  more  realistic  correlation  [5]. 
Since  the  wakes  of  compressor  blades  are  relatively  thick  compared  to  the 
trailing  edge,  the  difference  between  the  correlations  using  the  different 
parameters  is  considerable,  a  factor  of  four  in  frequency.  For  the  blade 
section  studied  here,  a  frequency  of  15  kHz  is  predicted  using  the  wake  thick¬ 
ness,  matching  both  the  experimental  and  numerical  results. 

The  modulation  of  the  frequency  and  strength  of  the  vortex  shedding 
observed  in  the  CFD  calculation  is  extremely  important  for  the  interpretation 
of  experimental  measurements.  The  presence  of  the  modulation  considerably 
complicates  the  practical  problem  of  vortex  detection  in  rotating  machinery. 

The  impact  on  laser  anemometry  techniques  and  spectral  analysis  methods  needs 
tc  be  quantitatively  assessed. 

Examination  of  the  details  of  the  numerical  simulation  as  well  as  classi¬ 
cal  vortex  shedding  analysis  [16]  suggests  that  the  process  may  be  considerably 
more  complex  than  can  be  represented  by  a  simple  Strouhal  number.  Many 
shedding  processes  may  compete  and  the  separation  zone  observed  in  the  simuia- 
tion  certainly  plays  a  part.  A  more  complete  discussion  of  vortex  formation 
and  evolution  will  be  the  subject  of  a  later  paper.  v. 

This  work  has  tended  to  treat  the  compressor  blade  and  vortex  street  in 
isolation  wner.  in  reality  it  is  part  of  a  very  complex  environment.  Acoustic 


feedback  has  been  shown  to  both  enhance  shedding  and  alter  the  frequency  [7], 
[13],  [19].  There  are  many  mechanisms  involving  inter-blade  row  interactions 
for  forcing  the  shedding  as  well.  The  possibility  of  phase  locking  between 
blades  is  also  evident  [17].  Clearly  much  work  can  be  done  in  this  area. 

Laser  anemometer  measurements  of  the  turbulent  kinetic  energy  within  the 
stator  passages  of  a  transonic  compressor  stage  have  shown  a  structure  consis¬ 
tent  with  a  vortex  street  wake  (i.e.,  periodic  islands  of  high  kinetic  energy 
in  the  wake)  [18].  A  laser  anemometer  measures  the  stationary  frame  ensemble 
averaged  flowfield.  Thus,  the  vortices  must  be  shed  at  fixed  phase  relative 
to  the  stator  in  order  to  appear  in  the  data.  This  implies  that  there  is  an 
additional  mechanism  by  which  the  measurement  of  time  averaged  performance 
using  closely  coupled  instrumentation  (in  the  stator  leading  edge,  for 
example)  may  give  misleading  results.  The  assumption  that  the  flowfield  is 
well  sampled  would  no  longer  be  valid.  The  stage  geometry  is  such  that  a 
stationary  frame  probe  would  sample  along  a  vertical  line  in  the  flow  maps 
given  in  Figs.  6-8.  Phase  locking  implies  that  the  patterns  evident  in  the 
figures  are  stationary  in  the  flow-wise  direction  with  respect  to  the  measure¬ 
ment  station.  Thus,  the  probe  will  not  be  uniformly  surveying  the  flowfield. 
The  quantitative  importance  of  this  effect  has  yet  to  be  determined. 

Another  area  needing  more  investigation  is  the  three-dimensional  wake 
structure  in  real  turbomachines.  Although  it  is  often  dangerous  to  generalize 
from  two  to  three  dimensions,  we  will  point  out  that  if  the  vortices  have 
considerable  spanwise  extent,  they  can  serve  as  an  extremely  powerful 
mechanism  for  the  radial  transport  of  f luid--complicating  the  evaluation  of 
the  radial  work  and  efficiency  distributions.  Work  with  cylinder  shedding  has 
shown  that  external  forcing  (acoustic  in  this  case)  can  enhance  the  spanwise 
coherence  of  the  shed  vortex  street  [19]. 
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A  change  ir.  measured  performance  due  to  the  wake  structure  is  a  result  of 
the  difference  between  common  probe  (generally  time)  averages  and  true  mass 
averages.  Thus,  in  addition  to  the  vortex  wake  structure  discussed  in  this 
paper,  other  processes  which  redistribute  thermal  energy  in  the  flow  (and  thus 
change  the  difference  between  the  averages)  can  have  similar  effects.  Local 
heat  transfer,  either  within  the  flow  or  to  the  blades,  may  be  important  in 
this  regard.  Furthermore,  the  dynamic  behavior  of  aerodynamic  probes  can 
introduce  artifacts  into  the  performance  measurements  which  could  either 
amplify  or  reduce  the  effects  noted  above,  depending  upon  the  exact  structure 
of  the  flow. 

In  this  paper  we  have  made  the  following  observations  concerning  the 
structure  of  compressor  blade  flow  and  wakes: 

1.  The  wakes  of  high  Reynolds  number,  transonic  compressor  blades  can  consist 
of  shed  vortex  streets.  This  is  confirmed  by  measurement,  modelling,  and 
numerical  simulation. 

2.  The  vortex  shedding  scales  with  wake  thickness. 

3.  The  shedding  frequency  and  strength  are  sensitive  to  the  environment. 

4.  The  vortex  street  can  depress  the  average  wake  temperature. 

5.  The  importance  of  the  wake  structure  increases  with  the  square  of  the 
freestream  Mach  number. 

6.  Wake  structure  can  influence  the  measurement  of  compressor  efficiency. 

7.  Changes  in  wake  structure  can  be  mistaken  for  cha-.ges  in  compressor 
performance . 

8.  The  wake  structure  car.  drive  the  blade  shock  system,  inducing  loss. 

9.  The  wake  structure  car.  De  influenced  by  other  fluid  dynamic  instabilities 
present  in  a  blade  passage  which  themselves  may  be  important  for  perfor¬ 
mance  and  structural  dynamic  reasons . 


1C.  Caution  is  required  wner.  interpreting  steady  state  solutions  of  inherent! 
unsteady  fiowfieids. 

Overall,  vortex  shedding  in  transonic  compressor  blade  wakes  can  have 
significant  influence  on  compressor  behavior. 
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Fig.  1:  Transonic  compressor  rotor  exit 
total  pressure  near  hub:  (*, 
instantaneous  measurement,  (B) 
an  individual  blade  passage  as 
seen  once  per  revolution,  and 
tC)  an  ensemble  average  of  that 
individual  passage,  from  [5] 


Fig.  2:  baser  anemometer  measurement  of 
the  outflow  from  a  transonic 
compressor  r'tor:  (A)  average 
velocity  profile,  (B)  histogram 
of  velocity  statistics  in  core 
flow  showing  Gaussian  profile 
characteristic  of  turbulence, 
and  1C)  velocity  statistics  in 
wake  showing  bimodal  distribu¬ 
tion  attributed  to  vortex  street, 
from  [<]  (below) 
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9:  Instantaneous  vector  plot 

(arrow  orientation  is  velocity 
direction,  length  is  velocity 
magnitude)  of  the  trailing  edge 
region  of  a  transonic  airfoil 
calculated  with  a  2-D  time 
accurate  Navier-stokes  code. 
Note  the  separated  region  on 
the  suction  surface  and  the 
vortices  shed  into  the  watte. 
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Fig.  10:  Static  pressure  history  at 

trailing  edge  of  a  transonic 
airfoil.  Note  the  periodicity 
due  to  vortex  shedding  and  the 
modulation  of  that  shedding  by 
a  lower  frequency  disturbance. 
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LI:  The  numeric*!  simulation  shows 
•  correlation  between  the 
passage  shock  motion  and  the 
separation  point  movement 
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Fig.  12:  a  comparison  of  the  shedding 

frequency  predicted  by  the  CFD 
simulation  with  those  inferred 
from  measurements 
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A  comparison  of  the  rotor  rela¬ 
tive  total  pressure  time  history 
measured  m  a  transonic 
compressor  with  that  predicted 
by  the  CFD  simulation 


TASK  II:  EXPERIMENTAL  AND  THEORETICAL  STUDY  OF  FLOWS 
IN  COMPRESSOR  HUB/CASING  TREATMENT 

(Investigators:  E.M.  Greitzer,  C.S.  Tan,  N.  Lee) 

APPROACH 

This  project  is  a  combined  analytical/experimental  investigation  of  the 
flow  phenomena  responsible  for  the  substantial  enhancement  of  compressor 
stability  found  using  casing/hub  treatment.  The  work  is  focussed  on  detailed 
examination  of  differences  in  endwall  flow  features  associated  with  an  endwali 
configuration  that  has  worked  extremely  well  in  improving  stable  flow  range. 
A-further  issue,  much  wider  in  scope,  is  the  general  question  of  just  what  is 
"compressor  stall"  and  what  are  fruitful  approaches  to  delaying  its  onset. 

To  address  this  point,  f irst-of -a-kind  experiments  with  controlled  hub 
slot  suction  and  injection,  and  basic  computations  of  the  flow  in  the  grooves 
have  been  carried  out. 

EXPERIMENTAL  WORK 
Background 

During  the  past  year,  most  of  the  effort  has  been  spent  on  the  experimen¬ 
tal  work.  As  described  previously,  an  experimental  facility  was  designed  with 
a  cantilevered  stator  row  and  a  rotating  hub.  Blading  parameters  were  chosen 
so  that  the  stator  (rather  than  the  rotor)  stalled  first,  stall  occurred  at 
the  stator  hub,  the  stall  was  wall  stall,  and  the  experiment  was  set  up  to 
simulate  the  flow  over  the  tip  of  a  rotor. 

The  detailed  measurements  taken  with  hub  slots  showed  that  the  suction  at 
the  rear  of  the  stator  passage  was  very  effective  in  suppressing  the  growth  of 
blockage  and  hence  rotating  S"ail  |  1  ]  .  However,  the  use  of  a  closed  slot 
meant  that  there  was  no  way  to  separate  ti  e  presence  of  the  suction  from  the 
presence  of  the  jet  ou*:  of  the  front  of  "he  grooves,  and  hence  to  conclusively 

from  the  rear  or  high  velocity  flow  injection 


identify  which  mechanism. 


suet  ion 


from  the  front  of  the  slot,  was  more  important.  For  this  reason,  an  experi¬ 
ment  was  set  up  to  provide,  separately,  injection  or  suction.  The  experiments, 
which  are  described  in  the  next  section,  used  a  slotted  hub  open  to  an  inner 
plenum  chamber,  through  which  a  controlled  amount  of  suction  or  injection  flow 
could  be  fed. 

Overall  Experiment  Description 

There  are  two  parts  to  the  present  experiment.  Part  one  is  an  examination 
of  the  effect  of  injection  or  removal  flow  rate  on  the  pressure  rise  charac¬ 
teristics  (speedlines).  This  has  almost  been  completed.  Part  two  will  be 
focussed  on  the  detailed  stator  flowfield  measurement,  also  with  injection,/ 
removal.  This  will  take  place  after  the  completion  of  part  one. 

A  total  of  four  hubs  were  machined  for  the  experiment.  All  were  60 
degree  axial  skewed  (Fig.  1),  but  the  axial  length  of  the  slots  and  their 
location  on  the  hub  varied  (see  Fig.  2).  The  slots  extended  all  the  way 
through  the  hubs,  i.e.,  were  open  to  the  plenum.  The  hubs  originally  machined 
are  designated  as: 

Full  1  inch  slot 

1/4  inch  open  at  rear  (denoted  as  trailing  edge  in  figure) 

1/2  inch  open  at  front  (denoted  as  leading  edge  in  figure) 

1/2  inch  open  at  rear 

In  addition.  The  full  1  inch  slot  was  subsequently  modified  to  yield 
another  siot,  open  1/4  inch  at  the  front.  The  full  1  inch  slot  one  had 
slots  covering  the  middle  90%  of  the  axial  chord  of  the  stator.  The  1/4  inch 
open  at  rear  had  slots,  1/4  inch  long  axially,  close  to  the  stator  trailing 
edge,  etc. 

The  main  features  of  this  part  of  the  experiment  are  summarized  as  foLiows 


•  Aii  speediir.es  were  taken  at  a  compressor  speed  of  2600  rpm,  corresponding 


to  a  rotor  tip  Mach  number  of  0.24  and  a  Reynolds  number  based  on  blade 


chord  at  the  stator  midspan  of  1.0*10^.  (Previous  tests  showed  that  there 
was  little  Reynolds  number  sensitivity  for  the  baseline  hub  treatment  con- 
f iguration . ) 

•  The  smooth  hub  and  the  1  inch  treated  hub  were  first  tested  without 

inj ect. ion/ removal .  After  that,  a  run  was  made  using  the  1"  treated  hub 
with  the  bottom  of  the  slots  sealed  to  eliminate  interactions  between  the 
slots  and  the  plenum. 

•  Different  hubs  were  then  tested  with  the  injection/removal  rate  varied. 

The  maximum  rate  was  7%  of  the  main  flow.  With  the  1/4  inch  builds  this 
led  to  too  low  a  plenum  pressure,  and  damage  to  the  seal,  so  these  builds 
were  tested  only  up  to  suction  rates  of  3.5%.  For  reference,  the  measured 
mass  flow  (in  and  out)  through  the  original  hub  treatment  was  approximately 
3.5%. 

•  In  the  tests,  the  flow  near  the  stator  hub  was  monitored  by  hotwire;  this 
proved  helpful  in  identifying  the  stall  point. 

•  Control  experiments  were  carried  out  to  verify  that  leakage  through  the 
seal  had  no  effects  on  the  speedlines. 

Compressor  Performance:  Speedlines 

Speedlines  were  taken  for  all  of  the  hubs  as  described  above  and  the 
dependence  of  stator  static-to-static  pressure  rise  on  a)  compressor  flow 
rate,  and  b)  injection/removal  rate  documented.  In  this  report,  static 
pressure  rise  (DP)  is  non-dimensionalized  by  0.5  pU^,  where  U  is  the  mean 
rotor  blade  speed,  unless  otherwise  stated.  The  compressor  flow  rate  is 
presented  as  the  flow  coefficient,  Cx/U,  measured  upstream  of  the  inlet  guide 
vanes.  Injection/removal  rate  is  shown  as  percentage  of  the  compressor  f’. ow 


The  figure  of  merit  used  here  to  assess  the  degree  of  effectiveness  of 


any  treatment  is  the  non-dimensional  peak  static  pressure  rise  attained, 
DP/0.5pU2.  Although  this  can  be  debated  with,  the  basis  for  use  of  this 
criterion  is  that  it  is  the  endwall  regions  which  give  rise  to  the  stall,  and 
what  we  are  really  assessing  is  the  abiity  of  the  endwall  regions  to  cope  with 
a  given  pressure  rise.  (Note  that  relying  on  flow  coefficient  at  the  stall 
point  is  somewhat  ambiguous  because  the  flow  coefficient  was  measured  upstream 
of  the  stator.  Had  it  been  measured  downstream,  those  speedlines  with  flow 
removal  would  have  shifted  to  the  left  (that  is,  towards  lower  C  /U)  and  those 
with  flow  injection  to  the  right.  Normalizing  by  |pC2,  where  C  is  the  inlet 
velocity,  also  carries  some  ambiguity  because  of  this.  For  these  reasons  we 
have  adopted  the  convention  just  described. ) 

Baseline  Runs 

The  baseline  runs,  which  are  the  smooth  wall  build  and  the  1  inch  treated 
build  (sealed),  are  discussed  first.  The  smooth  wall  speedline  is  displayed 
in  Fig.  3,  together  with  the  two  cases  tested  with  the  1  inch  treated  hub. 

The  smooth  wall  plot  is  similar  to  that  obtained  previously,  stalling  at 
Cx/U  =  0.337  and  DP/0.5pU2  =  0.083. 

The  same  figure  shows  results  from  the  1  inch  treated  hub,  sealed  to 
decouple  the  flow  in  the  slots  from  the  plenum.  For  this  case,  stall  occurred 
at  Cx/U  =  0.303,  with  peak  DP/0.5pU2  =  0.137. 

With  the  1  inch  treated  hub  unsealed,  the  stall  was  delayed  to  an  even 
lower  Cx/U  of  0.295,  a  13%  decrease  in  Cx/U  compared  to  the  smooth  wall  case. 
The  peak  DP/0.5pU2  attained,  as  well  as  the  abruptness  of  the  drop  at  stall, 
was  the  same  as  before. 

Speedlines:  No  Injection  or  Removal 

Figure  4  presents  the  speedlines  obtained  with  no  fluid  injection  or 


removal.  Compared  to  the  smooth  wall  build,  the  1/4  inch  hubs  gave  an 
improvement  of  5%,  whereas  the  1/2  inch  builds  achieved  36%  and  38%.  The  drop 
in  performance  of  the  1/2  inch  slots  at  stall  was  similar  to  that  of  the  1 
inch  slots. 

Speedlines  With  Flow  Removal 

Figures  5,  6  and  7  show  speedlines  of  the  1/2  inch  and  1/4  inch  builds 
with  different  suction  rates.  The  speedlines  are  arranged  such  that  those  for 
the  same  build  are  shown  together  in  one  of  these  three  figures.  The  baselines 
are  included  for  comparison. 

Compared  to  the  smooth  wall  characteristic,  it  appears  that  stall  was 
suppressed  regardless  of  where  (front  or  rear)  or  how  (1/2  inch  or  1/4  inch) 
this  suction  was  applied.  The  amount  of  increase  in  peak,  pressure  rise  was 
roughly  the  same  for  a  given  mass  extraction  independent  of  the  slot  geometry. 
With  suction  at  the  rear,  however,  performance  generally  dropped  abruptly  at 
stall.  This  was  not  so  when  suction  was  applied  at  the  front. 

Speedlines  With  Flow  Injection 

Speedlines  obtained  with  different  injection  rates  are  shown  in  Figs.  8 
through  10.  Injection  at  the  stator  hub  increased  the  stalling  pressure  rise, 
with  injection  at  the  front  of  the  passage  being  more  effective  than  injection 
at  the  rear.  Injection  at  the  rear  resulted  in  a  rather  abrupt  drop  in  per¬ 
formance  at  stall,  whereas,  with  injection  at  the  front,  the  effect  of  the 
stall  onset  on  the  compressor  characteristic  was  very  mild. 

Effect  of  Slot  Flow  Rate  on  Peak  Pressure  Rise 


In  Fig.  11,  stalling  pressure  rise  is  plotted  with  in jection/removal  rate 
for  all  of  the  hubs.  The  results  are  first  of  a  kind  which  shows  the  effect 
of  mass  flow  rate  and  injection/suction  location  on  stall  onset. 

For  the  1/4  inch  slots,  the  stalling  pressure  rise  increased  monotonically 
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with  suction  rate  in  the  whole  suction  domain.  For  the  1/2  inch  slot  tested. 


however,  this  was  not  the  case.  With  no  suction,  the  increase  in  stalling 
pressure  rise  of  the  latter  was  38%.  As  suction  increased  from  null,  there 
was  an  initial  drop  in  performance,  reaching  a  lowest  value  at  about  1.5% 
suction,  after  which  it  improved  with  no  further  drop.  The  zero-suction  value 
was  reattained  at  2.2%  suction. 

All  four  slots  had  closely  similar  performance  for  suction  higher  than 
2%,  with  the  performance  at  the  level  of  the  1  inch  treated  hub  (sealed)  at 
4.5%  suction. 

Using  DP/0.5pU^  as  a  criterion  with  blowing,  the  slots  at  the  front  of 
the  passage  performed  better  than  those  at  the  rear  at  the  same  injection  rate 
throughout  the  injection  domain.  For  the  1/2  inch  build,  the  "dip"  in 
performance  occurred  ir  uch  the  same  way  as  with  suction.  Note  that,  even  at 
6%  blowing  rate,  the  stalling  pressure  rise  was  lower  than  that  achieved  by 
the  1  inch  sealed  treated  hub. 

Work  During  the  Next  Contract  Period 

Although  there  are  some  points  to  fill  in  on  our  study  of  the  effect  of 
injection  and  removal,  this  part  of  the  experiment  can  be  regarded  as  almost 
complete.  The  next  major  task  is  to  take  the  detailed  hot  wire  measurements 
in  the  stator  passage  to  compare  with  the  hub  treatment  measurements.  This 
will  be  done  over  the  next  two  months.  In  addition,  we  plan  to  spend  consi¬ 
derable  time  in  analysis  of  both  the  hot  wire  and  the  overall  results. 

It  is  to  be  emphasized  that  all  of  the  results  presented  above  are  rela¬ 
tively  recent,  and  we  have  not  yet  analyzed  them  in  any  sort  of  depth.  Some 
preliminary  conclusions,  however,  can  be  drawn.  First,  suction  does  appear  to 
be  slightly  more  effective  than  injection  from  the  fron^  part  of  the  passage. 


and  considerably  more  effective  than  injection  in  the  rear  part  of  the  passage. 


for  flow  rates  up  to  that  measured  with  the  original  hub  treatment.  For 
higher  suction/injection  flow  rates,  the  suction  appears  to  have  a  signifi¬ 
cantly  stronger  effect  on  peak  pressure  rise. 

Second,  although  it  is  clear  that  the  suction  is  quite  important,  as 
hypothesized,  it  is  also  clear  that  the  actual  situation  is  much  more  complex 
than  can  be  explained  by  suction  alone,  and  that  with  the  treatment  slots, 
the  injection  near  the  front  end  is  also  of  use. 

Beyond  this,  there  are  other  issues  which  we  have  not  yet  addressed. 

Key  among  these  is  the  basic  question  of  what  phenomena  comprise  compressor 
stall.  In  particular,  the  hot  wire  results  from  the  previous  experiments 
with  the  hub  treatment  showed  that  the  fluid  on  the  hub  had  a  high  dynamic 
pressure  near  the  rear  of  the  stator  passage,  but  that  the  direction  of  the 
flow  was  basically  tangential,  with  a  very  low  axial  velocity.  Where  did 
this  fluid  come  from,  and  what  determines  the  velocity  field  near  the  wall? 
Questions  such  as  this  are  not  at  all  resolved,  but  are  at  the  heart  of  any 
type  of  rational  analysis  of  stall  and  flow  field  instability. 

Analytical  Work 

During  the  past  year  we  have  extended  the  simple  analysis,  reported  on 
previously,  for  the  flow  in  the  grooves.  The  overall  conclusions  were  similar 
to  those  reported  before  concerning  the  comparisons  of  computation  and  experi¬ 
mental  measurements  for  mass  flew  into  and  out  of  the  grooves,  and  hence  will 
not  be  discussed  in  any  detail.  The  main  point  is  that  the  analysis  is  now  on 
a  more  rigorous,  less  ad  hoc,  basis.  The  level  of  effort  on  this  work  has 
decreased  with  the  departure  of  Dr.  Modi,  although  a  detailed  report  on  the 
work  is  in  preparation. 

The  next  step  in  the  analysis  of  the  interaction  of  the  flow  in  the 


grooves  and  that  in  the  stator  (or  rotor)  passage  should,  we  think,  make  use 


of  the  three-dimensional  computational  procedures  that  are  now  being  developed 
Two  computations  in  particular  suggest  themself.  First  of  all,  a  blade 


passage  with  a  rotating  endwall,  such  as  is  in  the  present  experiment,  as  well 
as  in  rocor  tips.  Second,  some  modification  to  the  calculation  to  try  to 
simulate  the  inflow  and  outflow  from  the  grooves,  in  order  to  see  if  the 
observed  trends  in  peak  pressure  rise  versus  groove  geometry  can  indeed  be 
calculated.  We  hope  to  attempt  this  as  a  collaborative  effort  with  Allison 
personnel . 

Reference 

1.  Johnson,  M.C.  and  Greitzer,  E.M.,  "Effects  of  Slotted  Hub  and  Casing 


Treatments  on  Compressor  Endwall  Flow  Fields,"  to  be  published  in  ASME 
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TASK  III:  COMPUTATIONAL  TECHNIQUES  FOR  TURBOMACHINES 
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Progress  on  ISES  Design  System 


In  the  last  year  we  have  made  a  number  of  extensions  to  ISES.  Wher 


analysing  isolated  airfoils,  the  program  can  now  handle  supersonic  freestreams 


in  addition  to  subsonic.  Although  this  appears  to  be  a  minor  change,  it  has 


actually  been  technically  non-trivial  since  it  involves  the  question  of  what 


happens  to  the  Kutta  condition  when  the  trailing  edge  flow  becomes  supersonic. 


In  this  case,  there  is  no  physical  mechanism  for  ensuring  that  there  is  no 


pressure  jump  across  the  trailing  edge,  and  in  general  there  will  be  such  a 


jump  with  a  shock  and  expansion  fan  ensuring  that  there  is  no  pressure  jump 


across  the  wake.  In  normal  time-marching  codes,  there  is  numerical  viscosity 


which  in  some  ill-defined  sense  mimics  the  natural  viscosity  and  establishes 


the  "correct"  condition.  In  our  scheme  we  have  no  shear  viscosity,  but  have 


some  bulk  viscosity  in  supersonic  regions.  At  subsonic  speeds  we  have  to 


explicitly  impose  the  Kutta  condition  in  exactly  the  same  way  as  potential 


tetnods.  Surprisingly,  we  also  have  to  impose  the  Kutta  condition  at  super¬ 


sonic  speeds,  but  in  this  case  what  we  obtain  numerically  is  a  "boundary-layer" 


type  behavior  near  the  trailing  edge,  and  a  few  cells  away  from  the  trailing 


edge  there  is  an  appreciable  pressure  jump  across  the  airfoil.  When  analysing 


:ascades  the  situation  becomes  further  complicated  because  the  inflow  and 


outflow  can  be  in  any  one  of  three  states:  1)  subsonic,  2)  supersonic,  but 


axially  subsonic,  or  3)  axially  supersonic.  Each  different  state  requires 


different  boundary  treatments,  making  the  programming  implementation  much 


sarder  than  for  t  ime-jnarchir.g  methods,  but  perr.aps  more  rigorously  correct. 


In  addition  -o  the  -wo  papers  presented  at  the  AIAA  4th  Applied  Aerodyna¬ 


mic  Conference  in  San  Diego  in  .July,  which  nave  been  sen-  previously,  one 
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paper  wi.il  be  presented  at  the  AIAA  Aerospace  Sciences  Meeting  in  Reno  in 
January.  This  discusses  the  application  of  the  method  to  a  series  of  transonic 
viscous  airfoil  test  cases  which  has  been  collected  for  a  workshop,  organized 
by  Dr.  Pulliam  of  NASA  Ames,  which  is  being  held  as  part  of  the  meeting.  An 
advance  copy  of  the  paper  will  be  sent  as  soon  as  it  is  completed. 

Progress  on  the  research  monograph  has  been  slower  than  intended  due  to 
time  limitations.  Currently,  we  are  about  half-f inished  and  the  aim  now  is  to 
finish  over  the  Christmas/New  Year  period  when  there  are  fewer  other  commit¬ 
ments  . 

The  ISES  program  is  now  being  used  in  industry  by  Boeing  (commercial, 
military  and  helicopter  divisions).  General  Electric  (Lynn  and  Cincinnati  aero 
engine  divisions),  Lockheed  (Skunk  Works),  Northrop  ( R&D  center),  and  licensing 
discussions  are  being  held  with  United  Technologies  (Sikorsky  and  Pratt  & 
Whitney)  and  J.  Roncz  (an  independent  consultant  who  worked  with  Rutan  in 
designing  the  Beech  Starship).  The  program  is  also  being  used  by  government 
research  groups  at  AFWAL  (Dr.  A.  Wennerstrom) ,  NASA  Langley,  and  David  Taylor 
Naval  R&D  Center  (Dr.  T.  Huang). 

Progress  on  New  Research  on  Unsteady  Flow 

As  discussed  in  the  Research  Progress  and  Forecast  Report,  we  are  now 
redirecting  our  new  efforts  into  developing  methods  for  calculating  unsteady 
flows  in  turbomachinery  and  improving  our  understanding  of  these  flows  in 
conjunction  with  the  experimental  work  being  performed  by  the  other  investiga¬ 
tors  in  the  Gas  Turbine  Laboratory  and  their  associates  elsewhere.  At 
present,  ’-.he  computational  work  is  progressing  in  two  areas. 

The  first  area  is  the  development  of  a  new  algorithm  for  unsteady  tran¬ 
sonic  flows,  particularly  those  in  which  unsteady  shcck/boundary  layer 
interaction  is  one  of  the  Drin.cioal  features. 


Related  computational  work  in 


the  GTL  has  shown  that  the  vortex  shedding  in  transonic  compressors  is  due 
largely  to  a  separated  boundary  layer  on  the  suction  surface  which  in  turn  is 
caused  by  a  normal  shock  in  the  blade  passage.  Low  frequency  oscillations  of 
the  shock  have  been  calculated  by  an  earlier  method,  at  a  frequency  which  is 
dangerously  close  to  a  blade  bending  frequency  and  has  aroused  interest  in  the 
engine  industry.  It  is  believed  that  this  may  be  an  indication  of  a  shock/b.l 
instability  similar  to  those  found  in  transonic  diffusers. 

The  new  algorithm  is  an  extension  of  flux-vector-splitting  ideas  in  which 
the  gradient  of  the  flow  variables,  as  well  as  the  flow  variables  themselves, 
are  stored  at  the  cell  centers.  Using  the  gradients  allows  one  to  obtain 
second  order  accuracy  even  on  irregular  meshes,  while  using  the  flux-vector- 
splitting  allows  sharply  defined  shocks  without  requiring  excessive  grid 
resolution.  At  present,  the  algorithm  has  been  developed  for  inviscid  flows, 
but  in  the  future  it  will  be  extended  to  the  full  Navier-Stokes  equations  with 
the  intention  of  calculating  shock/b.l.  interactions.  A  paper  on  this  new 
method  is  being  submitted  to  the  AIAA  Computational  Fluid  Dynamics  Meeting  in 
Hawaii  next  June.  A  copy  of  the  extended  abstract  is  attached. 

The  second  area  is  the  calculation  of  turbomachinery  flows  with  vortices. 
Under  funding  from  Rolls-Royce  and  ONR,  another  program  has  been  developed 
which  uses  Ni's  Lax-Wendroff  algorithm  to  calculate  the  inviscid  flow  in  a 
stator  or  rotor  row.  The  boundary  conditions  in  this  program  allow  one  to 
specify  incoming  wakes  from  an  upstream  blade  row  and/or  incoming  isentropic 
pressure  disturbances  from  an  upstream  or  downstream  blade  row.  One  feature 
is  the  ability  to  analyze  cases  in  which  the  two  blade  rows  have  different 
pitches . 

The  Rolls  Royce/ONR  funding  is  primarily  for  algorithm  development,  and 
so  under  this  multi- investigator  grant,  in  parallel  to  the  experimental  work 


being  performed,  we  intend  to  use  this  program  to  investigate  the  importance 
of  vortex  shedding  on  the  performance  of  the  next  blade  row.  This  wili  be 
done  by  first  calculating  the  unsteady  flow  due  to  a  wake/rotor  interaction, 
and  then  repeating  the  calculation  but  instead  modelling  the  incoming  wake  as 
an  unsteady  Karman  vortex  street  instead  of  as  a  smooth  wake.  This  will  allow 
us  to  see  how  important  it  is  to  have  the  correct  wake  model,  which  is  clearly 
of  interest  to  both  the  theoretical  analysts  who  wish  to  use  simple  mathemati¬ 
cal  models  and  the  experimentalists  who  wish  to  assign  an  order  of  priority  to 
the  many  interesting  questions  to  be  addressed. 
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Abstract 
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Introduction 


Time-marching  schemes  for  the  Euler  equations  have  become  commonplace.  The  majority 
of  such  schemes  have  been  optimized  for  steady  calculations,  but  more  recently  unsteady  Euler 
simulations  have  received  increased  interest.  In  many  cases,  adequate  unsteady  simulations  cannot 
be  made  by  removing  the  time-inaccurate  mechanisms  (eg.  residual  smoothing,  multigrid,  etc.) 
from  a  steady-state  time-marching  scheme.  This  is  because  the  objectives  of  an  unsteady  simulation 
are  different  enough  to  require  a  new  formulation  of  the  solver. 

For  unsteady  Euler  simulations  it  is  necessary  to  have  correct  wave  propagation  prediction 
This  is  best  accomplished  by  schemes  which  attempt  to  directly  model  the  underlying  physics  of 
the  Euler  equations.  One  such  class  of  schemes  is  Riemann  solvers  [4],  but  these  are  difficult  to 
extend  to  general  2-D  grids.  Another  class  of  schemes  uses  flux-vector  splitting.  Included  in  this 
class  are  schemes  by  Steger  and  Warming  [5)  and  Anderson,  Thomas  and  van  Leer  [1],  These 
methods  use  the  correct  domain  of  dependence  for  flux  calculations,  which  is  more  important  in 
unsteady  flow  than  for  steady  calculations.  In  supersonic  flow,  central  difference  schemes,  such 
as  Jameson  et  al  [2],  and  partially  upwinded  schemes,  such  as  Ni  [3],  have  an  enlarged  domain  of 
dependence  and  hence  allow  information  to  propagate  upstream. 

The  Euler  equations  are  nondissipative,  hence  it  is  desirable  to  have  a  scheme  with  minimal 
dissipation.  This  is  accomplished  by  achieving  high  spatial  and  temporal  accuracy  (at  least  2nd 
order).  The  majority  of  time-marching  schemes  in  use  today  assume  a  Cartesian  grid  in  their 
formulation  and  are  2nd  order  accurate  on  such  grids.  However,  when  they  are  used  on  non- 
Cartesian  grids  their  order  of  accuracy  is  reduced,  resulting  in  increased  dissipation.  For  example, 
the  scheme  in  Ref.  [2]  uses  central  differencing  of  cell  averages.  Turkel  [6]  has  analytically  shown 
very  strict  restrictions  on  grid  smoothness  for  achieving  2nd  order  accuracy  with  this  scheme  in 
1-D.  He  also  showed  that  the  scheme  becomes  inconsistent  for  highly  stretched  grids.  The  schemes 
of  Refs.  [1]  and  [5]  use  upwind  differencing  of  cell  averages,  and  hence  also  suffer  from  these  order 
of  accuracy  problems. 

On  the  other  hand,  Ni’s  node-based  scheme  [3]  uses  trapezoidal  integration  to  calculate  the  flux 
balance  on  each  cell.  For  steady-state  calculations  this  is  spatially  2nd  order  accurate  (neglecting 
smoothing)  regardless  of  grid  smoothness.  Unfortunately,  spatial  accuracy  is  reduced  for  unsteady 
flows  because  the  distribution  step  (the  discretization  of  the  time  derivative  term  of  the  Euler 
equations)  assumes  a  Cartesian  grid. 

The  goals  of  the  present  investigation  are  two-fold.  First  to  develop  an  unsteady  time-marching 
scheme  which  gives  good  wave  propagation  prediction  through  the  use  of  flux-vector  splitting. 
Second  to  reduce  dissipation  by  achieving  grid  independent  2nd  order  accuracy  in  both  space  and 
time. 
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Technical  Approach 


2-D  Euler  Equations 


The  flow  of  an  inviscid  compressible  fluid  is  governed  by  the  Euler  equations.  They  represent 
conservation  of  mass,  momentum,  and  energy.  For  2-D  unsteady  flow  the  integral  form  is 


~  [[  UdA+cf  (F  dy  -  G  dx)  =  0, 

at  JJci  Jan 


(1) 


where  the  integrals  are  taken  over  a  control  volume  Cl  bounded  by  the  curve  dCl.  The  conservation 
variables  U  and  the  Cartesian  flux  functions  F  and  G  are  given  by 
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where  p  is  density,  q  =  ui  +  vf  is  the  velocity,  P  is  pressure,  e  is  toteil  internal  energy  per  unit 
volume  and  /i0  is  total  enthalpy.  Closure  of  the  system  is  obtained  by  the  equation  of  state  for  a 
perfect  gas, 


P  =  (  7-1) 
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where  7  is  the  specific  heat  ratio  (for  air  7  =  1.40).  The  total  enthalpy  relation  is  then  given  by 
ph0  =  e  +  P. 


Spatial  Discretization 


The  spatial  discretization  for  this  scheme  originally  came  from  a  Godunov-type  method  by  van 
Leer  [7]  for  the  1-D  wave  equation.  Within  each  cell  the  solution  is  assumed  linear. 


U(x,  y)  =  U  +  [x  -  xc)Ux  +  (y  -  yc)Ur 


(4) 


'where  U  is  the  solution  average  over  the  cell,  and  (Ux, Uy)  is  the  average  gradient  for  the  cell 
These  may  be  associated  with  the  cell  centroid  (xc,  yc)  defined  by 
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With  this  intracell  distribution,  the  solution  is  not  necessarily  continuous  at  cell  faces. 


Substituting  this  linearized  solution  into  the  Euler  equations  and  integrating  over  the  cell  area 
gives  an  equation  for  the  time  evolution  of  the  cell  averages. 


~  (All)  +  £  H(Um)  •  nAs  =  0,  H  =  Ff 4-  Gj 


(6) 


where  A  is  the  cell  area,  n  the  outward  normal  and  As  the  length  of  each  face  of  the  cell.  The 
fluxes  through  the  faces  H(Um)  are  evaluated  at  the  midpoints  of  the  cell  faces  as  shown  in  Fig  1 
If  Um  is  known  to  2nd  order  at  each  face  midpoint,  then  Eq.  6  is  2nd  order  accurate  regardless  of 
grid  smoothness. 

This  finite  volume  formulation  of  the  discrete  Euler  equations  is  conservative,  which  gives  the 
correct  shock  jump  relations  and  shock  velocity  for  unsteady  flows. 

Flux-vector  Splitting 

The  fluxes  through  the  faces  are  obtained  using  flux-vector  splitting.  In  1-D  the  flux  F  is 
separated  into  two  components:  a  forward  flux  (F+)  and  a  backward  flux  (F~).  At  a  given  point, 
F+  is  evaluated  using  solution  information  extrapolated  from  the  left,  and  F-  is  evaluated  using 
right  extrapolated  information.  This  is  illustrated  in  Fig.  2  where  the  flux  through  the  face  at 
i  -t-  1/2  is  given  by 

'<♦./. = *♦<»;;,„> + f-iu.w  m 

where  the  extrapolated  solutions  are  obtained  using  average  and  gradient  information  from  the  two 
adjacent  cells. 

U-+1/J  =  Ui  +  {Xt  +  l/2  -  U,>i/2  =  u*+l  +  (xi+ 1/2  -  *i+l)Uxi+l  (8) 

These  extrapolated  solutions  are  2nd  order  accurate  as  long  as  the  gradients  are  known  to  at  least 
1st  order. 

The  flux  splitting  used  is  that  developed  by  van  Leer  [8j.  In  supersonic  flow  the  split  fluxes  are 
given  by 

F*  =  F,  F-=0  M>+ 1 

F+  =  0,  F"  =  F  M  <  - 1 

In  subsonic  flow  (|M|  <  l)  they  are  given  by 


F±  =  A±c  { 
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The  2-D  extension  of  this  scheme  is  accomplished  by  rotating  the  fluxes  normal  to  the  face 
The  splitting  is  then  based  on  the  normal  Mach  number  defined  as 

-  (q/c)  ■  n,  M,  =  {q/c)-s  (11) 
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where  n  is  the  outward  normal  to  the  cell.  For  completeness,  the  tangential  Mach  number  is  alsc 
given.  In  the  face  normal  coordinates,  the  subsonic  split  fluxes  are  given  as 

(  ±f(l±M„)’  \ 
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where  the  second  and  third  entries  are  normal  and  tangential  momentum,  respectively. 


Gradient  Equations 

There  are  several  options  for  calculating  the  gradients  (Ux,  Uy)  on  each  cell.  One  method 
is  to  central  difference  the  cell  averages.  This  is  simple  but  has  the  wrong  domain  of  dependence 
in  supersonic  flow,  allowing  information  to  propagate  upstream.  This  counteracts  the  use  of  flux 
splitting  to  calculate  the  cell  averages.  Also  this  method  is  unconditionally  unstable. 

Another  option  is  to  update  the  gradients  by  minimizing  the  solution  jump  at  cell  faces.  For  a 
uniform  grid  in  1-D  this  works  out  to  be 

VxT'  =  £  {ffiV,  -  iaxtr^,)  -  (D”_,  +  1axU^_,)}  (13) 

For  steady-state  calculations,  this  is  a  simple  and  efficient  approach;  however,  it  gives  poor  wave 
propagation  predictions.  Like  central  differencing  it  allows  information  to  propagate  upstream  in 
supersonic  flow. 

A  third  option  is  to  evolve  the  gradients  in  time  along  with  the  cell  averages  by  taking  the 
moment  of  the  Euler  equations  as  suggested  in  Ref.  [7],  If  properly  discretized  this  method  gives 
the  same  domain  of  dependence  as  the  average  equations,  which  leads  to  proper  wave  propagation 
prediction. 


-  F  =  0 


The  derivation  of  the  gradient  equations  is  first  shown  for  1-D.  The  moment  equations  are 
obtained  by  first  multiplying  the  1-D  Euler  equations  by  x. 

(d\J  dF  1  d(xV)  d(xF)  „  .  . 

‘{isr+ai"0}  ^  V +  V  - F  = 0  (14> 

Next,  the  linearized  solution  is  substituted  into  Eq.  14  and  integrated  over  the  cell  (x  =  —Ax/2  to  -f  Ax 
When  this  is  done,  the  cell  average  in  the  time  derivative  term  vanishes,  giving 

d  Ax3  +  Ax/2 

T7  -TTT^  +*F  -  AxF(U)  =  0  (15) 

dr  !<!  -Ax/2 

In  this  equation  the  third  term  is  obtained  by  assuming  F(U)  is  linear  within  the  cell. 


F(U)  =  F(U)  +  z—  TJx  +  0(Ax2) 
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To  2nd  order  accuracy  the  second  term  is  given  by 


(xF),i.1y2  =  (x,  +  1/2  -  X.)  {F  +  (U,  +  (xt>1/2  -  I|)U^,)  +  F  (U,t1  +  (xl  +  1/2  ~  X,  + 1 )U.v, * l )  ] 

(17) 

Note  that  the  term  in  brackets  is  the  flux  through  the  face  evaluated  for  the  cell  average  equations. 
Hence,  the  rather  expensive  operation  of  calculating  fluxes  through  the  faces  need  not  be  repeated 
for  the  gradient  equations. 

Ref.  [7]  notes  that  calculating  the  gradients  by  the  moment  equation  results  in  a  linearized 
approximation  which  gives  the  integrated  least  square  error  from  the  exact  solution  within  each 
cell. 

In  two  dimensions  the  x  and  y  moment  equations  result  in  a  coupled  system  for  the  cell  gradients 
(Uxi  tJy).  The  full  semi-discrete  scheme  in  2-D  is  then  given  by 


dt 
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where  Azz,  Ayy  and  Azy  are  second  area  moments  about  the  cell  centroid. 

Avv  =  JJix  ~  zc)JdA,  Azz  =  J J  (y  -  ye)1  dA,  Azy  =  JJ (x  -  xc)(y  -  yc)  dA  (19) 

Note  that  if  the  gradients  are  neglected,  the  resulting  scheme  is  identical  to  the  1st  order 
upwinded  optior  in  the  scheme  of  Ref.  [  1] . 

Time  Integration 

The  semi-discrete  equations  are  integrated  in  time  using  an  explicit  3-stage  Runge-Kutta  scheme. 
Given  the  nonlinear  O.D.E. 
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the  solution  u(t)  is  evolved  in  time  by 


u°  =  u" 

1  =u°+  aAf/° 
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u2  =u°+l/2At/1 
u3  =  u°  +  At/2 
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(21) 


un+l  =  u 


This  scheme  is  2nd  order  accurate  in  time  independent  of  the  value  of  the  first  stage  integration 
constant  a.  Thus  a  is  a  free  parameter  which  can  be  used  to  maximize  the  time  step  At. 


ol 


Stability  Analysis 


The  stability  of  this  scheme  is  analyzed  using  the  linearized  1-D  Euler  equations  rather  than 
the  simple  scalar  wave  equation.  This  is  because  the  greatest  stability  restriction  occurs  in  subsonic 
flow,  where  the  Euler  equations  have  mixed  eigenvalues.  Also,  the  effect  of  the  split  flux  vectors 
is  not  easily  obtained  from  the  wave  equation.  These  matters  are  complicated  by  the  coupling 
between  the  average  and  gradient  equations. 


The  linearized  1-D  Euler  equations  are  given  by 


au  au 

at  +  Adx  ~0, 


(22) 


A  single  stage  of  the  discretized  form  of  the  average  and  gradient  equations  on  a  uniform  grid  is 
then  given  by 

TT*  =  CF°  -  {  [a+(i  -  T-1)  +  A-(T  -  1)]  u*-1 
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AU*  =  AU°  -  6a~  {  [A+(l  +  T~l)  +  A ~{T  +  1)  -  2A]  U*-1 
+  \  [a+(1  +  T~x)  -  A'(T  +  1)]  ATI*-1} 
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(24) 


where  AU  =  AiUx  and  T  is  the  translation  operator  (TU,  =  UJ  +  i).  In  addition  the  split  flux 
Jacobians  are  defined  as 

±  dF± 
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The  Fourier  mode  to  be  examined  is 
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(26) 


where  Q  is  a  6  x  1  constant  vector.  Substitution  of  this  into  Eqs.  23  and  24  results  in  the  matrix 
equation 

Q"+I  =  GQn  (27) 

where  G  is  the  amplification  matrix  whose  eigenvalues  must  have  magnitudes  less  that  one  for 
stability.  For  a  given  Mach  number,  a,  9  and  CFL  number  defined  by 

A  =  (u^c)^  (28) 


the  eigenvalues  of  G  are  calculated  by  EISPACK. 

Figure  3  shows  the  effect  of  the  first  stage  integration  constant  a  on  the  stability  boundary  for 
the  1**  order  scheme  (Eq.  23  only  with  no  gradients).  For  supersonic  flow  the  maximum  CFL  is 
A  =  2  at  a  ss  0.17,  but  the  stability  boundary  is  sensitive  to  the  choice  of  ct.  For  subsonic  flow  the 
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maximum  CFL  decreases  to  A  =  1.9  at  M  =  0.  As  a  comparison,  the  maximum  CFL  for  a  3-stage 
central  differenced  scheme  is  A  =  2  independent  of  Mach  number. 

The  stability  boundary  for  the  full  2nd  order  scheme  (Eqs.  23  and  24)  is  shown  in  Fig.  4.  For 
supersonic  flow,  the  maximum  CFL  number  is  A  =  0.59  at  a  «  0.17  and  decreases  as  M  —  0. 
Again,  the  maximum  CFL  i3  sensitive  to  the  choice  of  a.  A  linear  curve  fit  for  the  maximum  CFL 
is  then  given  by 

/  0.51  +  0.08M  M<  1  ,9Q, 

m&x  "  0.59  M  >  1  (  ’ 


Shock  Capturing 


As  shown  in  Ref.  [8],  the  1“  order  upwinded  scheme  (ie.  averages  only)  captures  shocks  with 
at  most  two  interior  cells  with  no  additional  dissipation.  Unfortunately,  this  is  not  also  true  for 
the  full  2nd  order  scheme.  In  the  immediate  vicinity  of  shocks  the  gradients  must  be  limited  or 
turned  off.  This  is  done  by  adding  an  exponential  decay  term  to  the  gradient  equations.  If  Rx  is 
the  residual  for  the  x  moment  equation  in  1-D,  then  the  modified  equation  is 


1  Ax3/ 12 


Ax3/ 12  1 6*p\, 


The  decay  term  contains  a  switch  to  turn  it  on  only  in  the  vicinity  of  shocks.  This  switch  is  the 
second  difference  of  average  density.  In  smooth  flow,  \62p\  ~  C(Az2)  and  the  source  term  is  2nd 
order.  In  the  vicinity  of  shocks,  where  \6*p\  ~  1,  the  source  term  gives  zeroth  order  accurate 
gradients  and  a  1*‘  order  solution. 


Results 

In  this  section  results  are  shown  for  transonic  and  supersonic  flow  over  a  bi-circular  arc  cascade, 
also  known  as  Ni’s  Bump.  These  are  standard  test  cases  for  2-D  Euler  solvers.  In  addition,  an 
order  of  accuracy  study  is  conducted  for  the  present  scheme. 

Transonic  Ni  Bump 

The  first  test  case  is  a  transonic  20%  bi-circular  arc  cascade  with  an  inlet  Mach  number  of 
0.675.  The  computational  grid,  shown  in  Fig-  5,  is  64  by  16  cells  and  is  algebraically  generated 
The  solution  on  this  grid  is  shown  in  Figs.  6,  7  and  8.  Figure  6  is  a  plot  of  Mach  number  distribution 
on  the  upper  and  lower  walls  of  the  channel.  As  the  plot  shows,  the  shock  on  the  lower  wall  is 
captured  with  only  one  interior  cell.  Note  also  the  absence  of  any  pre  or  post-shock  oscillations  or 
overshoots.  Figure  7  shows  Mach  number  contours  for  the  channel  interior;  throughout  the  domain 
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the  solution  is  smooth.  Typically,  for  Euler  solvers  oscillations  in  density  near  shocks  and  ofher 
strong  gradients  are  much  more  acute  than  for  Mach  number.  Figure  8  shows  that  the  density  is 
also  smooth  throughout  the  domain,  even  near  the  shock. 

Supersonic  Ni  Bump 

The  same  series  of  plots  is  shown  for  a  supersonic  4%  bi-circular  arc  cascade  with  inlet  Mach 
number  of  1.4.  The  grid  for  this  case  (Fig.  9)  is  56  x  16  cells.  It  has  the  same  distribution  as  that 
for  the  transonic  case,  except  that  cells  have  been  deleted  from  the  forward  section  where  the  flow 
should  be  uniform.  The  Mach  number  plots  in  Figs.  10  and  11  show  oblique  shocks  generated  at 
the  leading  and  trailing  edges  of  the  arc,  which  reflect  off  the  upper  wall.  Fig.  10  shows  that  the 
lower  wall  shocks  are  captured  with  one  interior  cell.  However,  by  the  time  the  leading  edge  shock 
impinges  on  the  upper  wall,  it  is  smeared  with  three  interior  cells.  This  is  believed  to  be  due  to 
cross-stream  dissipation  introduced  by  the  split  fluxes. 

Order  of  Accuracy  Study 

The  order  of  accuracy  of  the  present  scheme  is  verified  ueing  a  completely  subsonic  channel 
flow  (Mini«t  =  1/2).  The  arc  of  the  previous  cases  is  replaced  by  a  10%  sin2  x  bump,  which 
smoothly  blends  into  the  straight  portions  of  the  lower  wall.  This  eliminates  leading  and  trailing 
edge  stagnation  points  and  their  associated  large  truncation  errors.  Also  no  shock  losses  occur 
since  the  flow  is  subsonic.  Hence,  all  losses  for  this  case  are  purely  numerical.  In  this  study,  RMS 
of  total  pressure  loss  over  the  domain  is  U3ed  as  a  measure  of  these  numerical  losses.  It  is  defined 
as 


Pt/PToo\\2=^ 

Ei.jU  ~  PTij/PtooV 

1 

(31) 

Pt  =  P(1  + 

VM2)^r 

(32) 

Total  pressure  loss  has  been  found  in  the  past  to  be  a  good  indicator  of  accuracy  for  Euler  solvers. 

Figure  13  shows  the  finest  grid  (64  x  16)  used  in  this  study.  Coarser  grids  (32  x  8,  16  x  4  and 
8  x  2)  are  obtained  by  deleting  every  other  grid  line  from  the  next  finer  grid.  As  before  these  grids 
are  algebraically  generated.  Figure  14  shows  the  losses  for  both  the  1**  and  2nd  order  schemes  along 
with  Jameson’s  central  difference  scheme  (2]  on  these  grids.  The  slope  of  each  line  is  the  scheme's 
order  of  accuracy.  The  slope  for  the  present  scheme  using  averages  only  is  0.63,  indicating  that 
in  this  norm  the  scheme  is  not  even  1“  order  accurate.  The  full  scheme  using  both  averages  and 
gradients  (Eqs.  18)  is  truely  2nd  order,  having  a  slope  of  2.1.  For  comparison,  the  central  difference 
scheme  is  also  plotted.  Even  though  the  grids  are  fairly  smooth,  this  scheme  fails  to  obtain  2nd 
order  accuracy,  since  the  slope  of  its  line  is  only  1.4.  It  is  expected  that  the  central  difference 
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scheme  will  further  decrease  in  accuracy  as  the  grid  is  made  less  smooth,  while  the  present  scheme 
will  remain  unconditionally  2nd  order  accurate. 

For  solutions  on  the  64  x  16  grid,  the  present  2nd  order  scheme  has  approximately  10°$  the  total 
pressure  error  of  the  central  difference  scheme  and  about  1%  that  of  the  1st  order  scheme  using 
averages  only. 

Summary 

The  new  developments  of  the  present  investigation  to  be  discussed  in  full  in  the  final  paper 
include  the  following: 

A  new  time-marching  algorithm  for  the  unsteady  2-D  Euler  equations  is  presented.  The  scheme 
utilizes  flux-vector  splitting  for  shock  capturing  and  correct  prediction  of  wave  propagation.  Grid 
independent  2nd  order  accuracy  is  achieved  in  a  formulation  where  averages  and  gradients  are 
stored  tor  each  cell.  This  more  closely  models  the  nondissipative  nature  of  the  Euler  equations,  by 
reducing  truncation  error  on  non-uniform  grids. 

A  stability  analysis  of  the  linearized  1-D  Euler  equations  is  conducted  for  this  scheme.  Using 
explicit  3-stage  Runge-Kutta  time  integration,  the  maximum  CFL  number  for  the  spatially  2nd 
order  scheme  is  found  to  be  approximately  1/4  that  for  the  spatially  1“  order  scheme.  In  subsonic 
flow,  the  use  of  flux-vector  splitting  results  in  a  slightly  lower  CFL  restriction. 

Results  for  steady  transonic  and  supersonic  channel  flow  are  presented.  The  final  paper  will 
include  studies  of  shock  motion  for  unsteady  channel  flows. 

References 

'lj  Anderson,  VV.  K.,  Thomas,  J.  L.,  and  van  Leer,  B.,  “A  Comparison  of  Finite  Volume  Flux 
Vector  Splittings  for  the  Euler  Equations,”  ALAA  Paper  85-0122,  January  1985. 

2;  Jameson,  A.,  Schmidt,  W.,  and  Turkel,  E.,  “Numerical  Solutions  of  the  Euler  Equations  by 
Finite  Volume  Methods  Using  Runge-Kutta  Time-Stepping  Schemes,”  AIAA  Paper  81-1259. 
June  1981. 

3  Ni,  R.  H.,  “A  Multiple-Grid  Scheme  for  Solving  the  Euler  Equations,”  AIAA  Journal,  Vol.  20. 
No.  11,  November  1981,  pp.  1565-1571. 

4  Roe,  P.  L.,  J.  Computational  Phya.  43  p.  357,  1981. 

5  Steger,  J.  L.,  and  Warming,  R.  F.,  “Flux  Vector  Splitting  of  the  Inviscid  Gasdynamic  Equations 
with  Application  to  Finite  Difference  Methods,”  J.  Computational  Phya.  40,  pp.  263-293,  1981 


10 


6  Turkel,  E.,  "Accuracy  of  Schemes  for  the  Euler  Equations  with  Non-Uniform  Meshes."  AIAA 
Paper  86-0341,  January  1986. 

Jj  van  Leer,  B.,  “Towards  the  Ultimate  Conservative  Difference  Scheme.  IV.  A  New  Approach 
to  Numerical  Convection,”  J.  Computational  Phys.  23,  pp.  276-299,  1977. 

iS]  van  Leer,  B.,  “Flux-vector  Splitting  for  the  Euler  Equations,”  ICASE  Report  No.  82-30, 
September  1982. 


.p"Tjrv 


C>6 


Figure  1.  Typical  conservation  cell 
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Figure  2:  Flux-vector  splitting  in  1-D 
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Figure  3:  Effect  of  first  stage  integration  constant  on  stability  boundary  of  1*'  order  scheme 


l-D  Euler  Stability  Analysis:  2nd  order  scheme  (averages  a..d  gradients) 
Stability  Boundary 
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Figure  4:  Effect  of  first  stage  integration  constant  on  stability  boundary  of  2nd  order  scheme 
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Transonic  Ni  Bump:  M  =  0.675,  f  —  10%,  64  x  16  grid 
Computational  grid 
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Figure  5:  Computational  grid  for  transonic  Ni  bump  (Af  =  0.675,  r  =  10%) 


Transonic  Ni  Bump:  M  =  0.675,  r  =  10%,  64  x  16  grid 
Mach  Number 


Figure  6:  Mach  number  distribution  on  upper  and  lower  walls  for  transonic  Ni  bump 
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Superaonic  Ni  Bump:  M  =  1.4,  r  =  i%,  56  x  16  grid 
Computational  grid 
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Figure  9:  Computational  grid  for  supersonic  Ni  bump  (M  =  1.4,  r  =  4%) 


Superaonic  Ni  Bump:  M  -  1.4,  r  =  4%,  56  x  16  grid 
Mach  Number 
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Figure  10:  Mach  number  distribution  on  upper  and  lower  walls  for  supersonic  Ni  bump 
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Supersonic  Ni  Bump:  M  =  1.4,  r  —  4%,  56  x  16  grid 
Mach  Number 
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Figure  11:  Mach  number  contours  for  supersonic  Ni  bump  (AM  =  0.02) 


Supersonic  Ni  Bump:  M  =  1.4,  r  =  498,  56  x  16  grid 
Deneity 

1  51 

inc=  0.020 


Y  1  0 


0.5H 


0  0 

0.0  0.5 

Figure  12: 


1.0  1.5  2.0  2.5 

X 

Density  contours  for  supersonic  Ni  bump  ( Ap  =  0.02) 
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Subtonic  tin3  z  Bump:  M  =  1/2,  r  =  10%,  64  X  16  grid 
Computational  grid 
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Figure  13:  Computational  grid  for  accuracy  study  (Af  =  1/2,  r  =  10%) 


Order  of  Accuracy:  Subaonic  tin  *  Channel 
RMS  of  Total  Pretture  Lott 
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Figure  14:  Order  of  accuracy  study:  Scheme  and  grid  effect  on  RMS  total  pressure  error 


TASK  IV:  THEORETICAL  MODELLING  OF  STABILITY  AND  UNSTEADINESS 
IN  TRANSONIC  COMPRESSOR  FLOW  FIELDS 

(Investigators:  J.E.  McCune,  E.  Imperatori) 

Introduction  (Foreward) 

This  report  continues  the  theme  of  its  predecessor  (November  1985)  in 
discussing  recent  progress  in  understanding  unsteady  effects  in  turbomachines, 
especially  in  transonic  compressors.  As  reported  previously,  unsteady  vortex 
patterns,  similar  to  "vortex  streets",  can  form  behind  a  blade  row.  These  were 
shown  to  lead  to  "temperature  separation"  --  a  total  temperature  depression  in 
the  time-averaged  wake  behind  each  blade  --  providing  a  partial  explanation  of 
recent  experimental  observations  in  transonic  machines.  Entropy  spottiness 
may  also  occur. 

(This  Task's  schedule  was  delayed  somewhat  in  the  past  months  due  to  time 
needed  by  the  principal  investigator  to  recover  from  surgery.  The  delay  is  cur¬ 
rently  being  made  up,  and  we  expect  to  be  back  "on  schedule"  before  Christmas.) 

Imbedded  Flows  and  Convected  Patterns 

In  steady  or  quasi-steady  flow  past  airfoils,  the  phenomenon  of  separation 
and  reattachment  of  the  boundary  layer  --  leading  to  so-called  separation 
"bubbles"  --  is  relatively  familiar.  The  flow  within  such  bubbles,  featuring 
closed  or  recirculating  streamlines,  generally  has  depressed  total  pressure 
(and  total  pressure  gradients)  and  represents  a  dispersed  but  contained  region 
of  vorticity.  This  is  an  example  of  an  "imbedded"  flow.  The  bubble  is 
separated  from  the  free  stream  by  a  strong  shear  layer  (usually  thin)  surround¬ 
ing  it,  so  that  the  circulation  calculated  around  any  path  within,  or  partially 
within,  the  bubble  is  not  constrained  in  the  usual  way  by  Kelvin's  theorem. 

The  "core"  of  a  tornado  is  another  example  of  an  imbedded  flow  of  this  type. 

Generally,  for  strongly  unsteady  flows,  key  properties  array  themselves 
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on  streaklines,  not  streamlines.  Moreover,  it  is  often  difficult  to  perceive 
of,  or  even  to  properly  define,  "recirculating"  streaklines.  But  in  a  well- 
developed  vortex  "street",  certain  patterns  persist  as  they  are  convected 
along  in  the  wake.  It  is  then  often  possible  to  move  into  a  reference  frame, 
"convecting"  with  the  pattern,  in  which  the  pattern  appears  to  be  steady,  or 
only  slowly-varying.  In  that  case,  in  the  new  frame,  streaklines  and 
streamlines  again  become  virtually  coincident.  It  is  then  relatively  easy  to 
define,  as  in  the  case  of  the  "bubble"  above,  imbedded,  but  convected,  flow 
regions . 

In  a  convecting  vortex  street  pattern,  well-defined  regions  of  intense 
but  dispersed  vorticity,  with  corresponding  entropy  (or  total  pressure)  varia¬ 
tions  appear  and  can  be  treated  as  "imbedded"  flows,  much  as  suggested  above. 
The  importance  of  this  fact  lies  in  the  recirculating,  or  cyclic,  feature  of 
the  flow  viewed  in  the  convecting  frame.  Because  of  the  cyclic  feature,  some 
key  properties  of  the  convecting  vorticity  pattern  can  be  deduced.  In  the 
following,  we  illustrate  this  feature  for  plane  2D  flow. 

The  Cyclic  Constraint  -  Vortex  "Core"  Properties 

Our  main  interest  in  this  Task  is  the  development  of  concentrated,  convec 
ting  vortex  patterns  in  blade  wakes  at  transonic  Mach  Nos.  Compressibility 
and  sound  wave  interaction  can  be  important.  But,  in  order  to  illustrate  the 
cyclic  constraint  feature  simply,  we  limit  our  discussion  to  the  low  Mach  No. 
situation  for  the  purposes  of  this  report. 

In  the  convecting  frame,  after  the  formation  phase,  the  flow  is  nearly 
steady.  The  vorticity  evolves  in  this  phase  according  to  the  dissipative 


Helmholtz  equation.  At  low  Mach  No.  in  plane  2D  steady  flow,  this  reduces  to 


where  q  is  the  fluid  velocity  in  the  convecting  frame,  q  =  |q|,  v  is  the 
kinematic  viscosity  and  is  the  z-component  of  the  vorticity,  i.e.,  7lz  = 
(curl  q)*e_.  "4"  is  a  coordinate  along  a  streamline  (or  streakline  in  the 

convecting  frame). 

If  p  is  the  streamfunction  so  that,  with  q  =  (u,v) 
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Here,  "n"  is  a  coordinate  across  streamlines,  orthogonal  to  "n"  and  "g." 

make  up  the  so-called  "natural  coordinates"  often  used  in  characteristic  theory. 
We  can  write 

=  f(*,l)  (5) 


and  note 


aS  *  £*\  ‘  q£* 


where  f^  =  3 f / 8ip .  We  wish  to  find  f(ip,d),  or  a  good  approximation  to  it, 
because  if  we  can  do  that,  we  can  determine  the  flow  by  solving 


V2iO  =  -f  ( <0,  2.  ] 


to  determine  t|>(x,y). 

Outside  the  thin  shear  layer  surrounding  the  imbedded  core  of  vorticity, 
we  expect  f(4),l)  to  depend  only  weakly  on  "1";  in  fact,  3f/32  =  0(v)  from  Eq. 
(1)  unless  V2f  is  itself  very  large.  We  note,  accordingly. 


Thus,  for  any  finite  v,  no  matter  how  small,  the  bracket  in  Eq.  (10)  must 
vanish.  Equation  (10)  is  limited  to  recirculating  streamlines  not  in  the 
shear  layer. 

Using  the  fact  that  f  in  this  region  in  only  weakly  dependent  on  l,  the 
distance  along  the  streamline,  we  can  replace  f(ij;,8,)  by  F(i|j)  in  the  integrals, 
obtaining  approximately 

0  =  -FF*  |  ^  +  F"  j)  qdf, 
or  the  constraint,  for  each  <(i 

i  •  r  ( ijj )  h 

m)F  (no  =  f  («i()  (id 

where  "F"  is  the  circulation  around  the  particular  streamline,  and  "t"  is  the 


corresponding  "drift  time": 


reflects  the  intensity  of  the  internal  flow,  "b"  emerges  as  a 


While  "u)0" 

distortion  parameter,  reflecting  the  rearrangement  of  vorticity  within  and 
adjustment  of  the  shape  of  the  core,  or  bubble,  as  it  interacts  with  the 
external  flow,  or  free  stream.  Of  course,  more  parameters,  reflecting  more 
details  of  the  interaction,  can  be  included,  consistent  with  Eq.  (11). 

Extension  of  these  ideas  to  higher  Mach  numbers  introduces,  of  course, 
significant  complexity  due  to  compressibility  but  the  basic  usefulness  of  Eq. 
(11)  or  its  generalization  remains. 

Variations  from  Radial  Equilibrium 

Use  of  the  cyclic  constraint  above  and  the  determination  of  the  corres¬ 
ponding  core  flow  in  interaction  with  the  external  flow  (and  other  "cores") 
will  enable  us  to  evaluate  the  effectiveness  of  earlier  assumptions  concerning 
radial  equilibrium.  With  the  correctly  distorted  core,  an  improved  description 
of  the  momentum  balance  becomes  available  and  more  accurate  descriptions  of 
both  the  total  temperature  and  entropy  distributions  are  then  possible. 

Plans  for  Future  Work 

The  ideas  described  above  are  being  applied  to  the  complex  case  of  the 
wakes  in  a  transonic  compressor  environment.  3D  effects  are  to  be  included, 
using  techniques  indicated  in  earlier  reports. 

Reference 

1.  Batchelor,  G.K.,  "On  Steady  Laminar  Flow  With  Closed  Streamlines  at  Large 
Reynolds  Number,"  Journal  of  Fluid  Mechanics,  Voi.  1,  1956,  p.  177. 


V>v\V, 


3.  PUBLICATIONS  AND  PRESENTATIONS 

W.F.  Ng  and  A.H.  Epstein,  "A  Quasi-Three-Dimensional  Model  for  Intra- 
Stator  Transport  of  Rotor  Wakes,"  submitted  to  ASME  Journal  (presented  at  ASME 
Winter  Annual  Meeting,  November  1985). 

E.M.  Greitzer,  "Flow  Instabilities  in  Turbomachines,"  Review  Article  in 
Thermodynamics  and  Fluid  Mechanics  of  Turbomachinery,  Martinus  Nijhoff 
Publishers  (proceedings  of  NATO  Advanced  Study  Institute),  January  1986. 

H.W.  Shin,  E.M.  Greitzer,  W.K.  Cheng,  C.S.  Tan,  and  C.L.  Shippee,  "Circu¬ 
lation  Measurements  and  Vortical  Structure  in  an  Inlet  Vortex  Flow  Field," 
Journal  of  Fluid  Mechanics,  Vol.  162,  January  1986,  pp.  463-487. 

C.S.  Tan  and  E.M.  Greitzer,  "Compressible  Swirling  Flow  in  Turbomachine 
Annuli,"  AIAA  Journal,  Vol.  24,  No.  1,  January  1986,  pp.  92-101. 

H.W.  Shin,  W.K.  Cheng,  E.M.  Greitzer,  and  C.S.  Tan,  "Inlet  Vortex 
Formation  due  to  Ambient  Vorticity  Intensification,"  AIAA  Journal,  Vol.  24, 

No.  4,  April  1986,  p.  687. 

M.  Drela  and  M.  Giles,  "Viscous-Inviscid  Analysis  of  Transonic  and  Low 
Reynolds  Number  Airfoils,"  AIAA-86-1786-CP,  presented  at  the  AIAA  4th  Applied 
Aerodynamics  Conference,  June  9-11,  1986. 

M.  Giles  and  M.  Drela,  "A  Two-Dimensional  Transonic  Aerodynamic  Design 
Method,"  AIAA-86- 1793 -CP,  presented  at  the  AIAA  4th  Applied  Aerodynamics 
Conference,  June  9-11,  1986. 

A.H.  Epstein,  J.B.  Gertz,  P.R.  Owen  and  M.B.  Giles,  "Vortex  Shedding  in 
Compressor  Blade  Wakes,"  presented  at  AGARD/PEP  Specialist  Meeting,  Munich, 
Germany,  September  1986;  also  submitted  to  AIAA  J.  of  Propulsion  and  Power. 


M.  Drela  and  M.  Giles,  contributed  paper  to  Viscous  Transonic  Airfoil 
Workshop,  to  be  presented  at  the  AIAA  25th  Aerospace  Sciences  Meeting,  January 
12-15,  1987. 

M.C.  Johnson  and  E.M.  Greitzer,  "Effects  of  Slotted  Hub  and  Casing  Treat¬ 
ments  on  Compressor  Endwall  Flow  Fields,"  presented  at  1986  ASME  International 
Gas  Turbine  Conference,  Dusseldorf,  Germany,  June  1986;  also  to  be  published 
in  ASME  J.  Eng.  Gas  Turbines  and  Power. 

C.S.  Tan,  "A  Mixed  Spectral-Finite  Difference  Scheme  for  Three-Dimens  Iona 
Viscous  Flow  in  a  Bend  of  Rectangular  Cross-Section,"  submitted  to  Journal  of 
Computational  Physics. 

M.  Kurosaka,  J.  Gertz,  J.E.  Graham,  W.L.  Hankey,  and  P.  Sundown,  "Energy 
Separation  in  a  Vortex  Street,  Part  II:  The  Mechanism,"  to  be  published  :r. 
Journal  of  Fluid  Mechanics. 

M.  Giles  and  M.  Drela,  Two-Dimensional  Aerodynamic  Design  and  Analysis, 
to  be  published  by  Cambridge  University  Press  in  1987. 


4 .  PROGRAM  PERSONNEL 
Principal  Investigator: 

Edward  M.  Greitzer 

Professor  of  Aeronautics  and  Astronautics 
Director,  Gas  Turbine  Laboratory 

Co- Investigators: 

Alan  H.  Epstein 
Associate  Professor  of  Aeronautics  and  Astronautics 
Associate  Director,  Gas  Turbine  Laboratory 

Michael  B.  Giles 

Assistant  Professor  of  Aeronautics  and  Astronautics 
James  E.  McCune 

Professor  of  Aeronautics  and  Astronautics 


Choon  S.  Tan 
Research  Associate 


Graduate  Research  Assistants  (total  contract  period): 

9/83  -  1/85  Mark  Johnson  (AFRAPT  Student)*,  "Effects  of  Hub 

Treatment  on  Compressor  Endwall  Flow  Fields" 

(thesis  title) 

9/81  -  6/85  Thong  Dang**,  "A  Two-Dimensional  Design  Method  for 

Highly-Loaded  Blades  in  Turbomachines"  (thesis  title) 

9/83  -  Present  Petros  Kotidis*,  "Investigation  of  the  Radial 

Transport  in  a  Transonic  Compressor  Rotor  Stage" 
(thesis  title) 

9/83  -  10/85  Mark  Drela**,  "Two-Dimensional  Transonic  Aerodynamic 

Design  and  Analysis  Using  the  Euler  Equations" 

(thesis  title) 

9/83  -  6/85  Michael  Giles**,  "Newton  Solution  of  Steady  Two 

Dimensional  Transonic  Flow"  (thesis  title) 

9/84  -  9/86  Phil  Owen*,  "Computational  Simulation  of  Unsteady 

Flow  in  a  Transonic  Compressor  Rotor"  (thesis  title) 

9/84  -  Present  Norman  Lee 

9/85  -  Present  Edward  Imperatori 

9/86  -  Present  S.  Allmaras  (AFRAPT  student) 

*  S.M.  Degree  Completed 

**Ph.D  Thesis  Completed 


The  -as  TurD  me  ..aDoratory  a. so  r.as  an  active  seminar  program  to 
.•  rease  interact  ion  with  industry  and/or  government  by  bringing  speakers 
•  ■  MIT.  Turing  tie  period  covered  oy  'his  report,  these  have  included: 


Ar .  -'rarma.  Pra*  *  ■»  Wh  i  *  r.ey  Aircraft 
".mine  I,< is ses " 

W.  ;:'ee’<en.  >enera,  K  .ec*  r  !•'  'ompany 
’  -'  *  a  .  .  Phenomena  .  rs  Mu  .  *  is*  age  'ompressors" 

»:•  »  .  e.  .  e  r ,  ■•'a**  K  VI  h  :  *  r.ey  Ai rcraf * 

■  ««(•  -  Ihermogynam ; -:s” 

■  Vnt;pr  ,  '  .mm:  ns  Kr.g  :ne  '  .-mpany 

m.  .jprr  .  ,r,pr;  •  "  ,  rt)<  nar  ,er  ■  I  i  id  Mechanics" 

•  'r  w,  •>V  N  W  h  *  n  e  y  A.  rcraf*' 

.  it-*  .'is  Sew  -  •  ;.•.»»  -  •••  >r  H .  ar.K  :ece  >t  Paper 

■C  .  j ,r V  .  •  e ‘  .  •  :  ■  :  r" 


a  s  .  rr  r  «* 


■  i  ;»•  "  i  ■ 


jb„' 


7 .  CONCLUDING  REMARKS 


We  feel  that  the  multi-investigator  format  has  been  very  successful  in 
its  application  at  MIT.  One  facet  of  this  has  been  the  flexibility  to  assign 
resources  among  the  individual  projects.  Another,  perhaps  more  important, 
aspect  has  been  the  interactions  that  have  developed  between  various  investi¬ 
gators.  We  will  give  just  two  examples  of  such  interactions.  As  one.  Pro¬ 
fessors  Epstein  and  Greitzer  {plus  personnel  from  our  Guidance  and  Control 
Division)  are  starting  a  new  multi-disciplinary  research  area,  involving 
theoretical  and  experimental  investigation  of  suppression  of  turbomachinery 
instabilities  by  active  control.  As  a  second,  Phil  Owen,  a  student  of 
Epstein's,  was  able  to  use  Giles'  unsteady  Navier-Stokes  code  to  examine  a 
phenomenon  that  had  been  seen  experimentally.  This  latter  type  of  interaction 
where  state-of-the-art  computer  codes  are  used  to  derive  a  new  quantitative 
appreciation  for  phenomena  (rather  than  merely  run  as  design  tools)  seems  to 
us  to  be  at  the  core  of  this  type  of  multi-investigator  research  program.  A 
necessary  (although  not  sufficient)  condition  for  it  to  exist  is  the  existence 
of  a  critical  mass  of  faculty,  staff  and  students  who  are  interested  in  the 
type  of  complex  problems  that  characterize  high  performance  engines,  and  we 
hope  that  the  Gas  Turbine  Laboratory  can  continue  to  receive  support  to 
maintain  such  an  effort. 


